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ABSTRACT 


A  class  of  analytical  and  numerical  solutions  to  the 
problem  of  inviscid-jet  deflection  from  plane  and  axisym- 
metric  concave  surfaces  are  presented.   The  methods  of  solu- 
tion for  the  two-dimensional  cases  included  a  novel  method 
developed  in  the  present  investigation,  as  well  as  those 
due  to  Levi-Civita  and  Riemann  and  Hilbert.   These  methods 
have  provided  indirect  solutions  for  the  shape  of  the  curved 
solid  boundary  in  terms  of  a  given  jet-departure  and  jet- 
deflection  angle.   The  axisymmetric  cases,  namely,  the  deflec- 
tion of  axisymmetric  jets  from  hemispherical  thrust  reversers 
are  solved  directly  through  the  use  of  the  finite  element 
method  and  a  novel  iteration  scheme.   The  real-fluid  effects 
on  the  jet  deflection  are  estimated  by  comparing  the  results 
with  those  obtained  experimentally. 

The  analytical  and  numerical  methods  developed  or  used 
herein  are  sufficiently  general  to  yield  direct  and  indirect 
solutions  for  more  general  nozzle  and  curved-reverser  combi- 
nations of  special  interest  to  thrust  reversal  on  aircraft 
engines. 
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I.   INTRODUCTION 

The  prediction  of  the  deflection  of  finite  jets  by 
segmental  or  curved  deflectors  is  a  matter  of  importance 
in  the  design  of  target-type  thrust  reversers  on  jet-aircraft 
engines,  of  flip-buckets  on  spillways,  of  Pelton  wheels,  etc. 
In  general,  the  deflection  of  a  free  jet  by  a  solid  boundary 
is  well  suited  to  potential-flow  analysis  because  of  the 
dominance  of  inertia  and  pressure  intensity  in  the  establish- 
ment of  the  flow  pattern.   The  design  of  impulse  machinery, 
thrust  reversers,  etc.  utilizing  this  momentum  change  could 
be  facilitated  greatly  if  the  idealized  geometry  of  the  system 
under  potential  flow  conditions  were  known  because  such  condi- 
tions represent  asymptotic  values  which  are  approached  as  the 
effects  of  secondary  variables  such  as  entrainment,  boundary 
layer,  compressibility,  jet  attachment  to  adjacent  surfaces, 
etc.,  are  decreased.   With  such  information  available,  refine- 
ments of  design  -could  be  based  upon  a  secure  knowledge  of  the 
fundamentals,  and  many  rules  of  thumb  could  be  replaced  with 
precise  quantitative  data.   Specifically,  if  the  total  angle 
through  which  the  jet  is  deflected  is  determined  for  conditions 
of  both  partial  arid  complete  interception  by  the  boundary,  then 
the  principle  of  impulse  and  momentum  can  be  used  to  compute 
forces  or  other  dynamic  characteristics  of  the  system. 

The- problems  discussed  herein  fall  in  the  general  category 
of  finite  or  infinite,  two-dimensional  or  three-dimensional, 
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viscous  or  inviscid,  and  compressible  or  incompressible  flow 
past  bodies  of  arbitrary  shape.   In  fact,  in  its  broadest 
sense,  they  constitute  the  most  fundamental  and  also  most 
complex  problems  of  fluid  dynamics.   Cur  inability  to  deter- 
mine the  separation  points,  or  lines,  on  a  body  and  to  estab- 
lish  a  link  between  the  separation  and  the  unsteady  fluid 
motion  in  the  wake  of  the  body  gave  rise  to  numerous  exact 
and  approximate  methods  of  analysis.   Most  of  these  analyses 
rely  in  part  or  whole  on  the  experimentally  observed  charac- 
teristics of  flow  to  determine  the  magnitude  of  the  disposable 
parameters  introduced  into  the  analysis.   Suffice  it  to  say 
that  future  efforts  will  continue  to  be  based  upon  partly  a 
secure  knowledge  of  the  fundamentals  of  fluid  mechanics  and 
partly  on  experimental  facts. 

Evidently,  the  introduction  of  simplifying  assumptions  such 
as  inviscid,  incompressible  fluids,  the  use  of  bodies  with 
a  salient  edge  at  which  the  separation  point  is  fixed  and 
finally  the  restriction  of  the  analysis  to  either  two-dimen- 
sional or  axisymmetric  flows  enable  one  to  obtain  solutions 
which  are  often  of  practical  significance  and  in  conformity 
with  the  experimental  observations. 

The  present  work  deals  with  inviscid,  incompressible  flows 
of  finite  extent  impinging  upon  two-dimensional  or  axisymmetric 
bodies  with  sharp  edges,  and  makes  use  of  all  the  available 
exact  or  approximate  methods  of  analysis  in  arriving  at  novel 
solutions  towards  the  prediction  of  the  characteristics  of 
deflected  jets. 
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A  closer  look  at  both  the  problem  at  hand  and  the  methods 
available  requires  the  further  specification  of  the  geometri- 
cal configurations  of  both  the  nozzle  from  which  the  jet  is 
emerging  and  the  body  upon  which  the  jet  is  impinging.   The 
classical  theory  of  jet  deflection  deals  essentially  with 
nozzles  and  bodies  composed  of  straight  segments.   This  is 
primarily  because  the  successive  transformations  used  in  the 
analysis  lead  to  relatively  integrable  equations  for  the 
segmental  boundaries  and  to  exceedingly  complex  integral 
equations  in  the  case  of  curved  boundaries.   Furthermore , 
the  fact  that  there  cannot  be  a  corresponding  complex-function 
theory  for  three-dimensional  space,  the  limitations  of  our 
understanding  of  the  characteristics  of  axisymmetric  or,  in 
general,  three-dimensional  jet-deflection  problems  becomes 
quite  clear. 

In  order  to  bring  into  closer  focus  the  results  of  the 
past  efforts  as  well  as  the  limitations  of  the  existing 
methods,  we  will,  in  the  foregoing,  describe  briefly  the 
suggestions  made  concerning  the  "rounding  of  the  corners" 
of  solid  boundaries,  describe  the  essential  features  of  the 
problems  encountered  through  the  use  of  a  specific  example, 
and  then  discuss  the  history  of  the  two-dimensional  jet 
deflection  problem.   The  comments  concerning  the  previous 
work  on  axisymmetric  jet  deflection  will  be  reserved  for 
Section  III  where  such  cases  are  studied  through  the  use  of 
the  finite  element  method. 


A.   A  BRIEF  REVIEW  OF  SEGMENTAL  JET  DEFLECTORS  IN 
TWO-DIMENSIONAL  FLOWS 

The  two-dimensional  counterpart  of  the  jet-deflection 
problem  has  been  treated  by  several  investigators  through  the 
use  of  the  powerful  analytic-function  theory  and  successive 
conformal  transformations.   Sarpkaya  [1]  solved  the  U-shaped, 
two-segment,  deflector  problem  where  the  turning  angle  between 
the  segments  is  limited  to  90  degrees.   Tinney  et  al.  [2] 
extended  this  analysis  to  the  case  where  the  turning  angle 
between  the  symmetrically  situated  segments  is  greater  than 
90  degrees.   Later,  Chang  and  Conley  [3]  presented  an  analysis 
for  a  bucket  composed  of  a  series  of  segments  of  arbitrary 
number,  length,  and  angles;  however,  the  basic  as  well  as 
practical  problem  of  the  direct  analysis  of  jet  deflection 
by  curved  buckets  remains  unsolved  primarily  because  of  the 
mathematical  difficulties  encountered.   It  is  partly  for  the 
purpose  of  illustrating  these  difficulties  and  partly  for 
the  purpose  of  familiarizing  the  reader  with  the  fundamental 
concepts  and  the  analytical  transformations  that  the  basic 
steps  in  the  analysis  of  the  deflection  of  a  two-dimensional 
jet  from  a  U-shaped  segmental  deflector  will  be  discussed. 

The  flow  in  the  physical  z-plane  (see  Fig.  1)  may  be 
transformed  into  a  hodograph-plane  through  the  use  of 


„  _      ,       -i6   dW  ,,x 

C  =  -u  +  iv  =  -qe    =  -^  (1) 


The  flow  in  the  hodograph  plane  may  be  transformed  by  means 
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Figure  1   Transformation  planes  for  the  U-shaped  bucket 
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of  Planck's   logarithmic   transformation   into   an   ft-plane 
where    ft   is    defined  by, 


V. 
«   =   ln(-  -J-)    =   _   in    <L  +   ie  (2) 

6  J 


The  ft-plane  is  composed  of  straight  lines.   The  horizontal 
lines  represent  the  solid  boundaries  and  the  vertical  lines 
the  free  surfaces.   Such  a  polygon  may  be  transformed  (e.g. 
see  Milne-Tomson  [4])  through  the  use  of  the  Schwarz-Chris- 
toffel  transformation  onto  either  the  upper  or  the  lower 
half  of  a  t-plane  by  writing 


«   =    M0f 7^ 7-  +   N0  (3) 

°                          CU/7T                      dp/IT  ° 

(t-t1)      ±         (t-t2)      C     


Three  of  the  t  values  in  the  t-plane  may  be  chosen  arbitrarily 
in  accordance  with  the  Schwarz-Christoffel  transformation  and 
the  rest  must  be  assigned  unknown  parameters.   Accordingly, 
the  points  E,  A,  C  are  assigned  the  values  -k,  k  and  f 
respectively.   Then,  Eq.  (3)  becomes 


fi  =  Mo/  — ^      ?%     +   No  W 


(t2-k2)Vt^-l 


After  integration  and  evaluation  of  the  constants  M0 and  N^ 
one  obtains 


n  *  In  ^1  k  ,+  ,   3;"t —   for  t  >  1  and  t  <  -k 


y[^? 


(5) 

16 


and 


n  n  Vl-k2    +    V   1-t2  -  ,^ 

ft   =    In  for  -k   <_  t    <_  1  (6) 


The  potential  function  W  in  the  upper  half  of  the  t-plane 
may  be  written  in  terms  of  the  sources  and  sinks  representing 
the  flow.   Noting  that  there  is  a  source  at  D  and  a  sink  at 
C  and  that  only  the  upper  half  of  the  t-plane  is  used,  one 
has 


V.d 
W  =  -J—  [In  (t-f)  -  In  (t+1)]  (7) 

TT 


From  this  relationship  and  the  definition  of  r,   in  Eq.  (1) 

and  Q   in  Eq.  (2)  one  obtains  the  final  expression  that  relates 

the  z-plane  to  the  t-plane, 


Vl-t2  +  V  1-1  ' 


J  V  k2-t2 


Since  the  purpose  of  this  section  was  only  to  show  the  general 
procedure  used  for  this  type  of  transformations,  the  results 
of  the  integration  are  not  shown  here.   The  details  are  given 
by  Sarpkaya  [1]. 

It  is  apparent  from  the  foregoing  that  whenever  the  solid 
boundaries  in  the  physical  plane  (see  Fig.  1)  are  composed  of 
straight  segments,  then  the  ft-plane  also  is  composed  of 
straight  lines.   In  other  words,  the  flow  retains  a  constant 
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direction  along  the  solid  boundaries  (lines  parallel  to  the 
logarithmic  axis  in  the  fi-plane),  and  a  constant  velocity 
along  the  free  surfaces  (lines  normal  to  the  logarithmic 
axis  in  the  ft-plane).   Since  such  a  polygonal  boundary  can 
always  be  transformed  onto  a  t-plane  through  the  use  of  the 
Schwarz-Christof fel  transformation,  all  two-dimensional  jet- 
deflection  problems  of  this  nature  may  be,  at  least  theoreti- 
cally, solved. 

It  is  also  apparent  from  the  foregoing  that  whenever  both 
the  magnitude  and  the  direction  of  the  velocity  vary  along  a 
rigid  boundary,  then  the  fi-plane  is  not  in  general  composed 
of  straight  lines  or  of  lines  conveniently  intersecting  each 
other  at  90  or  180  degree  angle.   In  general  if  the  curved 
boundaries  in  the  physical  plane  is  given,  then  the  corres- 
ponding portion  of  the  ft-plane  will  be  an  unknown  curve 
since  the  magnitude  of  the  velocity  is  not  initially  known. 
One  may,  therefore,  attempt  to  simplify  the  problem,  for 
example,  by  specifying  the  variation  of  the  magnitude  of  the 
velocity  along  the  boundary.   This  does  not  in  any  way  alle- 
viate the  difficulties  to  be  encountered  later  since  there 
does  not  exist  a  transformation  function  comparable  to  that 
given  by  Schwarz  and  Christoffel  which  could  transform  an 
fi-plane  composed  in  part  by  curved  segments  onto  a  t-plane 
in  which  the  potential-function  may  be  written.   Consequently, 
one  will  either  seek  other  methods  of  handling  the  curved 
boundaries  or  completely  abandon  the  direct  approach  of 
obtaining  a  solution  for  a  given  geometry.   The  indirect 
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approach  called  into  action  will  seek  a  family  of  solutions 
for  a  family  of  initially  unknown  curved  boundaries  and  then 
will  let  the  designer  choose  one  among  those  obtained  to  suit 
his  desired  needs.   In  the  following,  the  efforts  made  by 
others  in  rounding  the  corners  of  buckets  otherwise  composed 
of  straight  lines  and  the  efforts  to  obtain  indirect  solutions 
will  be  described. 

B.   ROUNDING  OF  THE  CORNERS  AND  BOUNDARIES 

Because  of  its  theoretical  as  well  as  technological  im- 
portance, the  problem  of  rounding  the  sharp  corners  on  rigid 
boundaries  over  or  through  which  the  fluid  flows  has  attracted 
the  attention  of  many  mathematicians  and  hydrodynamicists. 
Prior  to  a  brief  exposition  of  these  efforts,  it  should  be 
noted  that  there  has  not  been,  up  to  now,  a  sufficiently 
general  method  of  studying  the  flow  around  curved  bodies  and 
that  the  techniques  so  far  developed  are  not  only  difficult 
to  handle  but  also  require  the  use  of  special  techniques  for 
a  given  problem.  • 

1.   Curved  Walls  of  Constant  Pressure 

As  cited  earlier,  both  the  magnitude  and  the  direc- 
tion of  velocity  vary  along  a  curved  boundary  and  the  func- 
tional relationship  between  the  wall  geometry  and  the  vector 
velocity  is  not  known  a  priori.   However,  by  making  a  simpli- 
fying assumption  such  that  the  pressure  or  the  total  velocity 
along  the  curved  wall  remains  constant,  one  can  round  the 
corners  in  a  bucket  otherwise  composed  of  straight  segments. 


19 


The  method  may  best  be  described  by  applying  it  to  the  U- 
shaped  bucket  case  discussed  earlier. 

Let  the  ft-plane  (see  Fig.  lc)  be  modified  in  such  a 
manner  so  that  it  takes  the  shape  shown  in  Fig.  2a.   Evidently 
the  magnitude  of  velocity  along  PQ  remains  constant  and  the 
direction  varies  from  tt/2  at  P  to  tt  at  Q.   The  ft-plane  is 
still  composed  of  straight  lines  and  may,  therefore,  be 
transformed  into  either  the  upper  or  the  lower  half  of  a 
t-plane  through  the  use  of  the  Schwarz-Christoffel  transfor- 
mation.  The  remainder  of  the  analysis  will  not  be  presented 
here.   Suffice  it  to  say  that  the  transformation  now  involves 
two  additional  terms  in  the  form  of  Wt-t   W  t-t   and, 
therefore,  renders  the  succeeding  integrations  relatively 
more  difficult  to  perform.   Be  that  as  it  may,  the  method 
is  in  principle  capable  of  generating  curved  corners  along 
which  the  pressure  or  the  velocity  remains  constant.   There 
is  sufficient  flexibility  in  the  analysis  to  move  the  points 
A  and  Q  along  the  boundaries  to  obtain  smaller  or  larger 
curved  corners  within  reasonable  limits. 

2.   Curved  Walls  with  Special  Pressure  Distributions 

A  bolder  approach  to  the  modification  of  the  ft-plane 
is  its  modification  in  such  a  manner  that  while  the  direction 
of  flow  along  the  curved  wall  varies  linearly,  the  magnitude 
of  the  velocity  varies  logarithmically.   This  procedure  still 
allows  the  ^-plane  to  be  composed  of  straight  segments  (see 
Fig.  2b).   It  is  immediately  apparent  that  the  application 
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Figure  2   Modifications  of  the  ft-  and  t-planes  to  round  the  corne 
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of  the  Schwarz-Christoffel  transformation  now  yields  terms 
such  as 


a  /tt       a  /tt 

(t-t  )  p    (t-t  )  q 
P  Q 


as  a  consequence  of  the  fact  that  the  line  segments  in  the 

ft-plane,  no  longer  intersect  each  other  normally.   Only 

under  rare  circumstances,  i.e.  for  special  values  of  a 

P 

or  a  ,  that  the  integral  equation  resulting  from  the  Schwarz- 
Christoffel  transformation  may  be  integrated  in  terms  of 
the  Beta- functions .   In  general,  the  transformation  may  be 
integrated  only  numerically.   Since  the  determination  of 
the  coordinates  of  the  rigid  boundaries  in  the  physical 
plane  requires  another  integration,  similar  to  Eq.  (8),  the 
solution  of  the  problem  via  a  two-stage  numerical  integration 
becomes  not  only  difficult  but  also  quite  approximate. 

This  method  has  been  applied  by  Schieldrop  [5]  to 
the  analysis  of  Borda-mouth-piece  type  flows  with  curved 
edges.   The  method  is  not  in  general  easy  to  apply  and  requires 
a  great  deal  of  ingenuity  and  sophistication  in  numerical 
integrations. 

3.   Curved  Walls  Through  the  Modification  of  the  t-plane 

Cockcroft  [6]  introduced  a  method  by  means  of  which 
the  corners  of  a  solid  boundary  may  be  slightly  rounded  by 
modifying  the  t-plane.   Consider  the  flow  in  a  90  degree 
corner  as  shown  in  Fig.  2c   This  plane  may  be  regarded  as 
either  z  or  the  fi-plane.   In  the  absence  of  round  corners, 
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the  Schwarz-Christof fel  transformation  relates  this  plane  to 
the  t-plane  (see  Fig.  2d)  by 


z  or  q  =  M„  f -ik 
-/  ft-; 


t-Q  dt 
(t-a)  Vt-b 


The  corner  at  B  may  be  rounded  in  an  indirect  way  by  replacing 
\J  t-0   in  the  above  equation  by  \j  t+1   +  Ay  t-1  .   In  this 
substitution  the  argument  of  V  t   which  changed  to  it/2  at 
t  =  0  now  varies  smoothly  when  t  moves  from  t  =  -1  to  t  =  +1. 
In  other  words,  a  curve  is  produced  in  the  z  or  fi-plane  whose 
curvature  depends  on  the  arbitrary  value  of  A.   It  should  be 
noted  that  for  a  value  of  A,  say 


b+1 
b-1 


the  velocity  or  the  pressure  along  the  curved  wall  remains 
constant.   In  fact,  evaluating  dW/dz  in  the  range  of  -l<t<l, 
one  obtains 


dW 


dz 


I  -I 

a 


This  method  was  successfully  used  by  Cockcroft  [6] 
in  calculating  the  electrical-stress  distribution  around 
certain  two-dimensional  conductors.   A  brief  discussion  of 
this  method  and  the  difficulties  associated  with  its  appli- 
cations are  also  described  by  Carrier  et  al.  [7]. 
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k .   Curving  of  the  Boundaries  Through  the  Use  of 
Special  Transformations 

The  Riemann-Hilbert  transformation  has  been  exten- 
sively used  by  Larock  and  Street  [8]  in  the  analysis  of 
flow  about  supercavitating  hydrofoils  and  by  Larock  [9]  In 
analyzing  in  an  indirect  manner  the  efflux  from  curved 
nozzles.   The  mathematical  details  of  this  method  will  be 
discussed  further  later  in  connection  with  its  application 
to  the  jet  deflection  problem. 

5 •   Analysis  of  Flow  about  Curved  Bodies  by 
Levi-Civita' s  Method 

Because  of  its  challenging  nature,  the  question  of 
the  existence  and  uniqueness  of  potential  flows  past  general 
obstacles,  having  boundaries  of  given  shape,  has  intrigued 
many  outstanding  mathematicians  such  as  Brillouin  [10], 
Villat  [11],  Levi-Civita  [12],  and  others.   In  fact,  Levi- 
Civita  may  be  said  to  have  solved  the  inverse  problem  of 
describing  the  class  of  all  jets  divided  by  curved  barriers. 
Although  no  specific  solutions  have  been  presented,  the  prob- 
lem of  determining  such  flows  was  reduced  in  special  cases 
to  the  solution  of  nonlinear  integral  equations  with  appro- 
priate boundary  conditions.   The  solutions  are  obtained  in 
general  by  solving  such  equations  numerically  and  by  suitably 
choosing  the  necessary  number  of  parameters.   Such  parameters 
described  for  example  the  variation  of  the  inclination  of 
the  velocity  vector,  the  variation  of  the  magnitude  of  veloc- 
ity, etc.   It  is  important  to  know  that  the  present  problem 
of  the  deflection  of  a  two-dimensional  jet  from  a  curved 
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boundary  has  not  been  previously  solved  by  any  of  the  methods 
so  far  described. 

In  Section  II,  Levi-Civita' s  method,  among  others, 
will  be  described  in  greater  detail  in  generating  a  family 
of  curved  buckets  for  a  given  jet  and  prescribed  jet 
deflection  angle. 

C.   SUMMARY 

Major  emphasis  in  the  foregoing  introduction  has  been 
on  the  direct  or  indirect  analysis  of  the  impingement  of 
ideal  fluids  on  two-dimensional  straight  or  curved  boundaries 
These  analyses  are  often  difficult  to  apply  and  the  amount 
of  work  required  may  not  be  commensurate  with  the  need  for 
a  solution.   In  such  cases,  the  use  of  relatively  less- 
exact  methods  may  be  adequate  or  even  necessary.   For  two- 
dimensional  flows  in  general  and  for  axisymmetric  flows  in 
particular  approximate  methods  have  been  developed.   A  de- 
tailed discussion  of  these  methods  with  the  exception  of  the 
finite  element  method  may  be  found  in  Robertson  [13].   These 
methods  may  be  classified  as  graphical  (field  plot  or  flow 
net),  finite-difference  numerical  analysis  (iteration,  relax- 
ation, Monte  Carlo),  and  analogies  (electrical  analogy,  Hele 
Shaw,  and  membrane  analogy),  and  finally  the  finite  element 
method.   The  first  three  of  these  approximate  methods  have 
been  widely  used  and  amply  described  by  Robertson  [13].   The 
finite  element  method  which  came  into  existence  during  the 
past  decade  is  currently  being  applied  to  the  solution  of 
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two-dimensional  and  axisymmetric  Laplace-field  problems. 

A  detailed  discussion  of  this  method  and  its  application  to 

axisymmetric  jet  deflection  problems  will  be  presented  in 

Section  III.   In  the  following,  exact  methods  of  analysis 

of  the  deflection  of  two-dimensional  jets  from  curved  buckets 

will  be  undertaken. 
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II.   ANALYSIS  OF  TWO-DIMENSIONAL  JET 
DEFLECTION  FROM  CURVED  BUCKETS 

In  the  following,  three  methods  will  be  employed  to 
obtain  a  family  of  curved-bucket  shapes  in  terms  of  a  given 
jet  deflection  angle  and  the  angle  of  departure  at  the  lip  of 
the  two-dimensional  bucket.   The  jet  will  be  assumed  to  be 
free  and  impinging  symmetrically  upon  the  bucket.   These 
methods  consist  of  the  modified  hodograph-method,  Levi- 
Civita's  method,  and  the  Riemann-Hilbert  method.   The  first 
of  these  is  novel  and  the  other  two  have  been  previously 
employed  as  discussed  in  Section  I.B. 

A.   MODIFIED  HODOGRAPH  METHOD 

The  method  consists  of  the  definition  of  the  hodograph- 
plane,  its  transformation  into  an  infinite  strip,  and  subse- 
quently, into  a  t-plane  through  the  use  of  the  Schwarz-Chris- 
toffel  transformation.   The  flow  in  the  t-plane  is  expressed 
in  terms  of  appropriate  singularities  and  then  the  geometry 
of  the  curved  bucket  is  determined  through  the  evaluation 
of  the  proper  integrals.. 

Consider  the  hodograph-plane  in  Fig.  3b  where 


r  =  -l^=±  e-10  (10) 

g     V,  dz   V.  K      J 

i  3 


The  fluid  stagnates  at  A,  accelerates  along  ANB  in  a  manner 
specified  by  the  hodograph-plane,  leaves  the  bucket  at  B  at 
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Figure  3   Transformation  planes  for  the  modif ied-hodograph  method 
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an  angle  of  71,  and  eventually  acquires  an  angle  of  deflection 
3.   Thus,  the  outer  unit  circle  in  the  hodograph-plane  repre- 
sents the  free  streamlines  of  constant  velocity.   As  stated 
earlier,  the  fluid  is  assumed  to  be  inviscid  and  the  real- 
fluid  effects  such  as  boundary  layer,  entrainment,  etc.  are 
ignored . 


Clearly,  the  arc  ANB  in  the  hodograph  plane  may  be  chosen 
to  represent  different  velocity  variations  along  the  bucket 
and  different  angles  of  departure  at  B,  e.g.,  the  curve  AE. 
Although  the  analysis  will  be  relatively  more  complex,  there 
are,  to  be  sure,  no  conceptual  difficulties  and  differences 
between  that  presented  herein  and  other  cases. 

The  hodograph-plane  may  be  transformed  to  an  S-plane 
(see  Fig.  3c)  through 


s  - .-  f-H  (11) 


and  the  interior  of  the  polygon  in  the  S-plane  may  be  trans- 
formed into  the  lower  half  of  the  t-plane  by 


S  =  -  In  t  (12) 

TT 


which  is  obtained  through  the  use  of  the  Schwarz-Christof fel 
transformation. 

The  flow  in  the  t-plane  is  comprised  of  two  sinks  and  a 
source.   Hence  the  complex  potential  is  given  by 
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V.d  , 

W  =  -J— [ln(t+k)  +  ln(t+r) 

it  k 


-  2  ln(t+l)]       (13) 


Writing 


dz  =  -  ~  -  dW 


to  obtain  the  geometry  of  the  bucket,  one  has 


dz  =  Sm  +  211^]  .  [1  +  1.2 

•ft  In   t  t+k        .  .1        t+1 

z  k 


]    dt  (14) 


After  separation  of  the  real  and  imaginary  part,  one  obtains 


Y   ,     (t+k)  (t+h 

§  =  -  In  ^~ 

d   *       (t+1)2 


(15) 


and 


L   = 
d 


-2  / 


(4-  + 


In  t  't+k     1   t+1 


)  dt 


(16) 


The  deflection  angle  R  is  related  to  the  parameter  k  through 


k  =  EXP  (-  tt  sln  3  ) 
1  -  cosR 


(17) 


— i  R 
which  is  obtained  by  writing  £  =  e    ,  inserting  it  in 

Eq.  (11)  and  combining  with  Eq.  (12). 
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The  foregoing  completes  the  formal  analysis  of  the  problem, 
The  bucket  geometry  may  now  be  calculated  for  a  given  deflec- 
tion angle  0  through  the  use  of  Eqs .  (15),  (16)  and  (17). 
The  results  obtained  through  a  straightforward  integration 
scheme  are  shown  in  Fig.  4  for  a  set  of  representative  jet 
deflection  angles.   The  coordinates  of  the  50  degree  case 
are  shown  in  Appendix  A,  and  the  computer  program  in  Appendix 
B.   Suffice  it  to  say  that  the  thrust  acting  on  the  bucket 
may  easily  be  determined  through  the  use  of  the  momentum 
equation  in  terms  of  the  deflection  angle  3. 

B.   LEVI- CI VITA'S  METHOD 

This  method,  which  is  fully  described  in  several  reference 
texts  (e.g.  Birkhoff  and  Zarantonello  [14],  Milne-Thomson 
[4]),  consists  of  the  definition  of  a  complex  function  co  as 


a)  =  6  +  i  In  (=^-) 


where  9  and  q  are'  respectively  the  direction  and  the  magnitude 
of  the  velocity  vector.   V.  is  the  constant  velocity  along 
the  free  streamlines.   The  co-function  is  expressed  in  terms 
of  a  polynomial  in  such  a  manner  that  w  is  real  over  the  free 
streamlines  and  complex  over  the  curved  obstacle.   The  shape 
or  the  curvature  of  the  obstacle  is  defined  in  terms  of  the 
coefficients  of  a  polynomial  representing  co. 

Even  though  this  method  was  originally  devised  to  deal 
with  curved  boundaries,  as  discussed  in  Section  I.B-5,  its 
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application  has  been  mostly  restricted  to  boundaries  with 
straight-line  segments  (Cisotti  [15]),  because  of  the  ana- 
lytical difficulties  encountered  in  its  application  to  curved 
boundaries.   Wu  [16]  has  employed  it  in  the  study  of  flow 
about  fully-cavitating,  slightly-curved,  hydrofoils. 

The  physical  z-plane  (see  Fig.  5a)  consists  of  a  jet  of 
width  2d  and  velocity  V.  impinging  upon  a  curved,  symmetric, 


two-dimensional  bucket  BAB1 .   The  complL^-potential  function 
W  is  given  as  usual  by  (see  Fig.  5b) 

W  *  <}>  +  iiji  (15) 

A  straightforward  application  of  the  Schwarz-Christof fel 
transformation  establishes  a  relation  between  the  t-  and 
W-p lanes .   Thus,  one  has 


V.d 

W  = ^—  [In  (1-t)  +  In  (1+t)]  (19) 

if 


Because  of  the  symmetry,  t  is  taken  equal  to  +1  and  -1  at 
C  and  C  respectively  and  zero  at  A.   The  parameter  k  is 
assigned  to  B  and,  as  will  be  shown  later,  it  is  uniquely 
related  to  the  deflection  angle  $. 

The  upper  half  of  the  t-plane  is  transformed  into  the 
inner  region  of  a  semicircle  of  unit  radius  (Fig.  5c)  through 
the  transformation  given  by  Levi-Civita  [12]  as 


t  =  -  |  (c  +  I)  (20) 
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Figure  5   Transformation  planes  for  Levi -Civita ' s  method 
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The  circular  arc  BAB'  in  Fig.  5c  represents  the  solid  boun- 
dary over  which  the  direction  of  the  velocity  vector  is 


assumed  to  be  prescribed  and  diametric  line  BCDC'B'  repre- 
sents the  free  streamlines  where  the  magnitude  of  the  velocity 
is  known . 

The  relationship  between  the  points  in  the  z,  t  and 
C-planes  are  given  by: 

point  B'  ,   t  =  -k,  X>   »  +1 

point  A  ,   t  =  0,    5  ■  i 

(21) 
point  D  ,   t  -*■  °°,    C  =  0 

point  B  ,   t  =  k,    c  =  -1 
The  W-plane  is  related  to  the  z-plane,  as  before,  by 


V4   dz    V.  e  K       J 


where  q  and  0  represent  respectively  the  magnitude  and  direc- 
tion of  the  velocity  vector. 

One  can  now  define  an  w-function  which  is  quite  similar 
to  that  given  by  Planck  as 


u  =  i  In  (-Jk  |)  .  e  +  i  Hi  (^L)  (23) 


Now  solving  for  z  in  the  above  equation  and  expressing  all 
the  variables  as  functions  of  c  one  has 
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d     V  ;        dr,  ^  K      ' 


Defining 


r,  -  r  e 


and  noting  that  over  the  bucket  only  the  argument  and  along 
the  free  streamlines  only  the  modulus  varies,  the  combination 
of  Eqs.  (19)  and  (24)  yields, 


2 
z  _  1   -  ioj(a)r2k   cosa  sina-,  ,  /ot-\ 

---./e     [- -g — ]  da  (25) 

1  -  k   cos  a 


and 


«  ■  I  .  /elu(r)[2;(r  "  r2) ]  dr  (26) 

"It  Ir  .     1  2 


Equation  (25)  gives  the  coordinates  of  the  bucket  and  Eq .  (26) 
the  coordinates  of  the  free  streamlines. 

The  function  oo  is  assumed,  according  to  Levi-Civita,  to 
be  of  the  form 


o)(c)  =  i  In  \l\   +  aQ  +  a^  +  a2c2  +  a^3  (27) 


where  the  logarithmic  term  can  be  identified  as  that  corres- 
ponding to  a  free  jet  impinging  normally  upon  a  flat  plate 
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(Cisotti[17l) •  The  additional  terms  in  the  polynomial  whose 
coefficients  are  to  be  prescribed  later  modifies  the  expres- 
sion for  cj  and  helps  to  curve  the  solid  boundary. 

Because  of  the  symmetry  of  the  flow  in  the  physical 
plane,  the  coefficients  a  ,  a„ ,  a^.,  ....  are  not  acceptable. 

Noting  that  over  the  solid  boundary  r  =  1  and  over  the 
free  streamlines  a  =  0  or  a  =  tt,  Eq.  (27)  may  be  reduced 
over  the  bucket  to 

/  n    .  n   I  cos  o\       L    it  ,     io  .     i3c      i5c 
w(a)  =  i  ln  1  +  sine  *  2  +  ale   +  a3e    +  V 

(28) 

where  plus  sign  is  to  be  used  for  tt  >  a  >  tt/2,  and  minus 
sign  for  tt/2  >  a  >  0.   Over  the  free  streamlines,  it  reduces 
to 


u)(r)  =  -  tan  -1  — ^—k   +  a,  r  +  a0r3  +  a^r5  .  .  .      (29) 

1  -  r2    1  3      5 


The  coefficients  a,,  a^,  a^,  ...  may  now  be  chosen  to  obtain 
different  bucket  shapes.   In  the  present  study,  a)  is  chosen 
to  satisfy  the  following  conditions 


CO 


(  +  1)   =  -TT 


co(i)   -  ±  p"  -  i  oo 

(30) 

0)(-l)   =  TT 


CO 


(o)   =  0 
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These  conditions  simply  state  that  the  jet  departure  angle 
at  the  lip  of  the  bucket  be  equal  to  ±ir,  and  that  the  bucket 
be  continuous  in  curvature  at  the  stagnation  point  A  and 
symmetric  with  respect  to  the  jet  axis. 

There  are  several  ways  to  select  the  coefficients  a. . 
In  this  study  a  Fourier-series  representation  was  used 
primarily  because  of  the  versatility  afforded  by  it  in 
representing  different  bucket  shapes. 


Equation  (28)  can  be  written  over  the  circular  arc  AB1 
as 


w(a)  =  i{ln  tj — r r^—   +  a,  sine  +  a_  sin  3a  +  a.  sin  5a...} 

i  t  sina    i  j  o 

(3D 
+{-  j  +  a,  cos  a  +  a~  cos  3a  +  a^  cos  5a  ...} 


The  real  part  of  Eq .  (31)  represents  the  direction  6  of  the 

velocity  vector  over  the  bucket.   Thus,  it  is  obvious  that, 

one  can  prescribe  the  curvature  of  the  bucket  or  the  modulus 

of  the  velocity  vector  by  properly  choosing  the  coefficients 

a.,  through  a  . 
1  n 

Let  6  be  represented  by  the  lines  shown  in  Fig.  6.   Evi- 
dently, by  choosing  sufficient  number  of  terms,  one  can  re- 
present the  variation  of  G  with  a  shown  in  Fig.  6.   In  doing 
so  it  is  important  to  make  0  symmetric  with  respect  to  the 
6  axis  since  the  series  are  expressed  in  cosine  terms  only. 

Rewriting  Eq.  (31)  one  has  over  the  arc  AB' 
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Figure  6  Approximate  representation  of  6(a) 
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w(a)  -  i  In  =  C,°l°  n   +   i  h(a)  +  0(a)        (32) 
i  t  sin  a 


where 


N 
h(a)  =   E   a  sin  na  (33) 

n=l  n 


and 


N 

6(a)  =  -  §  +   E   a  cos  na  (3*0 

d     t   n 
n=l 


where 


a  = ~  [cos  no,  -  cos  na0]         (35) 

1  ~  °2 

in  which  a-,  and  a?  represent  the  two  parameters  shown  in 

Fig.  6.   The  evaluation  of  the  Fourier  coefficients  a   is 
to  n 

shown  in  Appendix  C.  To  complete  the  analysis  one  needs  to 
evaluate  the  term  e  Q  which  appears  in  Eq .  (25)  in  terms 
of  Eq.  (32) .   Writing 


^iw(a)    1  +  sin  a  -h(a)   18(a) 

e      =  e      e 

cos  a 


and  separating  the  real  and  imaginary  parts,  one  has 


iw(a)    1  +  sin  a  -h(a)  /    rt/  N  ,  .   .  n/    w    r^c\ 

e   v  '  = e      (cos  6(a)  +  i  sin  0(a))  (36) 

cos  a 


kO 


Now  combining  Eqs.  (25),  (33),  (35)  and  (36),  one  finally 
obtains 


x   1     f  -h(a)   .   0/  x  2k   sin  a(l  +  sin  a)  , 

-r  =  —    /  e      sin  6(a)  ^ ~ da 

d   tt   /0  t    ,  2     2 

tt/2  1  -  k   cos  a 


(37) 


and 


y  _  1     ,   -h(a)     a,    v  2k  sin  a(l  +  sin  a)  . 

^-  =  —    /  e      cos  0(a)  s-* ?: da 

d   tt   /0  ...  2    2 

tt/2  1  -  k   cos  a 

(38) 

Equations  (37)  and  (38)  yield  the  coordinates  of  each 
segment  of  the  curved  boundary  through  the  use  of  the  appro- 
priate values  of  a-,    and  a~.   This  in  essence  completes  the 
application  of  the  Levi-Civita' s  method  to  the  analysis  of 
the  jet  deflection  from  curved  obstacles.   In  the  following 
the  implications  of  the  variations  of  a.,  and  a?  will  be 
discussed. 

The  variation  of  9  with  a  shown  in  Fig.  6  may  be  written 
for  each  interval  as 


0(a)    =   -tt  0   <   a   <   o-. 


6(a)    =   f   ( )    -  |   (2   +   ^ )  a±   <   a    <   a 

0(a)    =   -  f         c2<a<f  (39) 

IT  TT 

0(a)    =-p"         -^  <   a    <    tt  -  a~ 
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6(a)    =  f   (a2  -  a        "  a     -  S    >  "   -  °2   <   0    <    '  "  °i 


8(a)    =0  Tr-a<a<iT 


It  will  be  shown  that  this  particular  functional  depen- 
dence produces  bucket  geometries  of  the  form  shown  in  Fig.  7- 

Rewriting  Eqs .  (37)  and  (38)  for  the  particular  example 
under  consideration  one  has 


^  =  0  (40a) 


and 


i-  =  |    /  2  e-h(a)  k2  sln  "t1  +  |ln  0)  do   (HOb) 
tt/2  1  -  k  cos  a 


for  the  interval  of  r  >  o  >o?i    i.e.  for  the  segment  AP, 
(see  Fig.  7) 


*  =  2   /e-h(a)  k2  sin  a(l  +  sin  a)  cos[ir  °  "  ql  -,  da 

d   it  ,    .  ^    d         2  a0  -  a, 

a?  1  -  k  cos  a  2    1 

(4la) 

and 


J  =  2   ;Vh(a)  k?  slnp(^  +  s*n  q)  sin  [f  a  "  0l  ]  dr  +  yn 
d   *  a2  1  -  k2  cos2a         2  a2  "  al        p 

(41b) 
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for  the  interval  a   >  a  >  a  ,  i.e.  for  the  curved  segment 
PQ,  and  finally 


x  .  £   ;°  e-h(a)  k2  sin  (1  +  sin  a)  da  +         (i 
d   *  ax  1  -  k2  cos2a  q 

and 


5  =  y  (42b) 

d   Jq 


for  the  interval  a,  >  a  >  0,  i.e.  for  the  segment  QB. 

The  coordinates  of  the  free  streamlines  may  be  obtained 
by  integrating  Eq .  (26)  through  the  use  of  Eq.  (29)  and  the 
coefficients  a.  given  by  Eq.  (37).   Finally,  the  deflection 

angle  is  evaluated  by  replacing  in  Eq .  (29)  the  corresponding 

1        / 2~ 

value  of  r  at  the  point  C,  i.e.  r  =  r-(l-\l-k   ). 

The  foregoing  equations  have  been  numerically  integrated 
(the  computer  program  is  given  in  Appendix  D)  ,  and  the  re- 
sulting family  of  bucket  shapes  are  shown  in  Fig.  8a.   A 
special  case  for  a     =  0.31416,  o     =  0.47124,  and  k  =  0.946 
is  shown  in  Fig.  8b. 

C.   RIEMANN-HILBERT  METHOD 

The  method  consists  of  the  transformation  of  the  physical 
plane  into  the  W-plane,  where  W  is  given  by 


W  =  $  +  i  \Jj 
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and  the  transformation  of  the  V/-plane  into  the  upper  half 
of  a  t-plane  through  the  use  of  the  Schwarz-Christof fel 
transformation.   Then  the  flow  in  the  t-plane,  which  is 
composed  of  segments  of  solid  boundaries  and  free  stream- 
lines are  related  to  the  Riemann-Hilbert  transformation  as  a 
well-posed  mixed  boundary  value  problem.   Then  a  certain 
shape  function  f(t)  is  'assumed  and  the  characteristics  of 
the  flow  in  the  physical  plane  are  evaluated  through 
integration. 

The  details  of  the  technique  will  become  clear  as  it 
is  applied  to  the  determination  of  the  deflection  angle  of 
a  two-dimensional  jet  impinging  symmetrically  upon  a 
two-dimensional  bucket. 

Figure  9a  shows  a  jet  of  width  2d  impinging  upon  a  curved 
bucket  with  a  180  degree  angle  of  departure  at  the  upper  and 
lower  lip  of  the  bucket. 

It  is  apparent  from  the  application  of  Bernoulli's  equa- 
tion that  the  velocity  is  constant  along  the  free  surfaces. 
The  complex- function  W,  for  the  upper  half  of  the  physical 
plane,  is  represented  by  an  infinite  slit  as  shown  in  Fig.  9b 
The  W-plane  can  in  turn  be  transformed  into  the  upper  half 
of  a  t-plane  through  the  use  of  the  Schwarz-Christof fel 
transformation.   The  straightforward  application  of  this 
transformation  yields 


V,d 
W(t)  =  -i-  in  t  (42) 
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(a)    z-plane 


-**.  $ 


rmrrri  1 1 1  /  /  l 


(b)    W-plane 


0 
D 


1  k 

gzzzzzczzxczE, 


A 


(c)    t-plane 


Figure  9   Transformation  planes  for  the  Riemann-Hilbert  method 
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where  V.  is  the  velocity  along  the  free  surfaces.   In 
accordance  with  the  rule  of  this  transformation,  three  of 
the  t  values  are  assigned  arbitrary  values  (t  =  1  at  A, 
t  =  +°°  at  C  and  t  =  0  at  D) ,  and  an  unknown  parameter  k  is 
assigned  to  the  point  B. 

The  complex-potential-function  W  is  related  to  the 
physical  z-plane  by 


dW       -10  Mlos 


where  q  and  0  represent  respectively  the  magnitude  and  the 
direction  of  the  velocity  vector. 

It  is  convenient  at  this  point  to  separate  the  magnitude 
and  the  direction  of  the  velocity  through  the  application 
of  Planck's  logarithmic  transformation  given  by 


m»lnf  -ifl  (44) 

J 


Combining  Eqs.  (42)  and  (44)  and  solving  for  z  in  terms  of 
t,  one  has 


z  =  _  1  /  p~w(t)  dt  sue.) 

d      tt  ;  e       t         •  K^> 


It  is  apparent  from  the  foregoing  development  that  the  real 

part  of  the  complex-function  w(t)  is  zero  along  the  free 

streamlines  where  q  =  V..   Thus, 

J 
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Re[ui(t)]  =  0    for   t  >  k   and   t  <  0 

The  imaginary  part  of  w(t)  which  represents  the  direction 
of  the  velocity  vector  is  zero  along  the  line  AD  where 
0  <  t  <  1.   Over  the  bucket,  however,  both  the  magnitude 
and  the  direction  of  the  velocity  are  variable  and  are  an 
unknown  function  of  t".   Let  this  function  be  denoted  by  f(t) 
in  the  interval  1  <  t  <  k. 

The  foregoing  conditions  define  a  well  posed  Riemann- 
Hilbert  mixed-boundary  value  problem  where  the  real  and 
imaginary  parts  of  an  analytic  function  co(t)  are  alternately 
known  along  the  real  axis,  i.e., 

Re(co)  =  0    t  >  k 

Im(oj)  =  -f(t)  1  <  t  <  k  (46) 

Im(co)  =  0    0  <  t  <  1 

Re  (to)  =  0    t  <  0 

The  general  solution  of  the  Riemann-Hilbert  problem  is  well 
known  (Cheng  and  Rott  [18]  and  Larock  [9]). 
In  principle  one  writes 

w(t)  =  H(t)  Q(t)  (47) 
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where  H(t)  is  the  solution  of  the  homogeneous-Riemann  prob- 
lem (i.e.  for  f(t)  =  0)  as  it  satisfies  all  the  other  condi- 
tions listed  in  Eq.  (46)  and  may  be  chosen  as 


H(t)  =  ^-t(t  -  k)  (48) 

In  passing  it  should  be  noted  that  the  particular  form  of 
H(t)  which  constitutes  a  solution  for  f(t)  =  0  is  not  of 
special  importance  as  shown  by  Song  [19]  and  that  uniqueness 
of  the  solution  does  not  depend  on  the  particular  value  of 
H(t). 

The  function  Q(t)  may  be  written  as  (Larock  [91) 


Q(t)  =  I  ;+m   Im[Q(?)]  dn  (M9) 

7T       n  —  t 


Combining  Eqs.  (46),  (47)  and  (48),  one  has 


j  +  y        V-t(t  -  k)    rk    -f(n)       dn 

WltJ  - J        —  j— r-\ 

1  J   -n(n  -  k)  (n  "  t} 


(50) 


The  remainder  of  the  analysis  depends  on  the  selection 
of  f(t).   As  pointed  out  earlier,  f(t)  would  have  been  known 
had  the  shape  of  the  solid  boundary  been  known  a  priori.   The 
present  analysis,  however,  is  an  indirect  approach  and  seeks 
the  shape  of  the  boundary  in  terms  of  the  boundary  conditions 
assigned  and  the  deflection  angle  assumed.   Thus,  it  is 
apparent  that  one  may  select  f(t)  as  an  arbitrary  function 
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of  t  provided  that  the  boundary  conditions  at  the  stagnation 
point  and  at  the  lip  of  the  bucket  are  satisfied,  and  obtain 
a  family  of  bucket  shapes.   In  the  following,  f(t)  will  be 
taken  as  a  linear  function  of  t  partly  for  the  sake  of  sim- 
plicity and  partly  in  order  not  to  obscure  the  application 
of  the  method  by  introducing  rather  complex  and  unmanageable 
integrals.   Thus,  writing, 

f(t)  =  a  +  bt  (51) 

and  choosing  a  and  b  in  such  a  manner  that  the  slope  of  the 
bucket  be  ±  it/2  at  the  point  A  and  ±  it  at  the  lip,  one  has 


it  k  -  2 

a  =  7? 


2  k  -  1 


and 


b  -5  1 


•  2  k  -  1 


Inserting  Eq.  (51)  in  Eq.  (50)  one  obtains 


w(t)  =  V-t(t  -  k)   ;k  -(a  +  bn)   _dn_^     (52) 

1   V^ulT^lo  n  "  t 


The  integration  of  Eq.  (52),  the  details  of  which  are 
presented  in  Appendix  E,  yields, 
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Id 


(t) 


.  (a  +  bt)   -1 
i  cos 


-ft-t- 1 


t  -  1 


]-i|^0   COB-1^] 


(53) 


for  t  <  0,  and  t  >  k 


and 


u(t) 


-  <a  ;bt)  in  ' k(_t) ^(k-tlcos1^] 


(5*0 


for  0  <  t  <  k 

The  deflection  angle  3  can  be  found  by  letting  t  -*-  -°°  in 
Eq.  (53),  or  by  letting  t  -*■   °°  in  Eq .  (50)  and  then  performing 
the  integration  as  shown  in  Appendix  F.   The  deflection  angle 
3  is  found  to  be  uniquely  related  to  the  parameter  k  as 


3  =  it  + 


2(k  -  1) 


/3k  m   -12- 

(  4p-  2)C0S     — £■ 


k 


V 


k  -1 


(55) 


To  obtain  the  shape  of  the  bucket  and  ,the  actual  location  of 
the  free  streamlines,  Eqs .  (53)  and  (5*0  are  replaced  in 
Eq.  (45)  and  integrated. 

The  integration  must  be  performed  numerically  using  the 
method  described  by  Hartree  [20]  because  of  the  singularity 
at  the  point  t  =  1.   It  must  be  pointed  out  that  the 
singularity  at  t  =  1  is  an  integrable  singularity  of  the  form 
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/  1  $X 

o      v^r 


No  attempt  is  made  here  to  perform  the  numerical  integration 
of  Eq.  (45)  even  though  the  procedure  is  v,Tell  known  as  shown 
in  one  application  by  Tinney  [2],  because  the  method  previously 
described,  namely,  Levi-Civita' s  method,  is  very  similar  to 
the  one  presented  herein  and  has  the  advantage  of  being  very 
flexible  to  obtain  different  bucket  shapes  and  presents  no 
singularities  for  the  integration. 

D.   A  CRITIQUE  OF  THE  ANALYTICAL  METHODS  FOR 
TWO-DIMENSIONAL  FLOWS 

The  past  as  well  as  the  present  analyses  have  emphatically 

demonstrated  that  almost  all  of  the  difficulties  of  the 

analysis  of  jet  deflection  from  curved  obstacles  are  ascribable 

to  the  determination  of  the  initially  unknown  free  streamlines. 

In  view  of  this  fact  a  direct  analysis  of  the  problem,  that  is, 

the  determination  of  the  jet  deflection  and  the  position  of 

the  free  streamlines  via  either  purely  analytical  or  partly 

analytical  and  partly  numerical  techniques  is  at  best  an 

exceedingly  complex  proposition.   As  propounded  by  Birckhoff 

et  al .  [14]  and  other  hydrodynamicists  what  one  can  hope  for 

is  an  indirect  solution  of  the  problem;  that  is,  obtaining  a 

family  of  obstacle  shapes  which  will  yield  the  prescribed  jet 

deflection  angle  and  provide  a  familiarity  between  the  shapes 

so  obtained  and  the  range  and  change  of  direction  of  the  free 

parameters  involved  in  the  analysis.   What  has  been  said  so 
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far  holds  true  for  both  the  plane  and  axisymmetric  flows. 
Should  one  adopt  a  purely  numerical  procedure  (finite  differ- 
ence, relaxation,  finite  element,  etc.)  it  is  then  and  only 
then  that  it  is  possible  to  begin  with  a  prescribed  obstacle 
shape  and  to  approach  a  unique  solution  through  successive 
iterations . 

In  what  follows,  attention  will  be  focused  only  upon  the 
evaluation  of  the  individual  as  well  as  relative  merits  of 
the  analytical  methods  dealing  with  the  deflection  of  plane 
jets  . 

The  hodograph-method  developed  herein  is  capable,  within 
the  limitations  of  the  potential-flow  theory,  of  analyzing 
the  flow  in  curved  obstacles  and  generating  a  family  of 
practically  suitable  bucket  shapes.   Furthermore,  this  method 
has  the  inherent  capability  of  being  extended  to  the  genera- 
tion of  even  a  larger  family  of  obstacle-shapes  through  the 
use  of  distributed  sources  and  sinks  in  the  hodograph-plane . 
Even  though  this  concept  has  not  been  pursued  further  herein, 
there  are  no  conceptual  difficulties  in  doing  so. 

The  method  devised  by  Levi-Civita  and  successfully  em- 
ployed herein  requires  the  judicious  selection  of  the  varia- 
tion of  6  with  a  and  the  use  of  a  suitable  series  to  enable 
one  to  perform  the  integrals  leading  to  the  coordinates  of 
the  obstacle  as  well  as  of  the  free  surfaces.   The  deflection 
angle  may  be  determined  either  by  calculating  the  asymptotic 
slope  of  the  free  streamlines  or  by  finding  the  pressure 
distribution  over  the  bucket  and  using  the  equation  of 
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momentum.   The  emerging  fact  is  that  a  brief  familiarity  with 
this  method  enables  one  not  only  to  generate  a  family  of 
suitable  curved  obstacles  but  also  to  round  the  corners  of 
buckets  composed  of  straight  segments.   For  example,  by 
letting  (a.  -  a?)  -*■  0   one  finds  that  the  bucket  shape 
generated  by  this  method  exactly  approaches  that  previously 
solved  by  Sarpkaya  [1]  (see  Figure.  8a).   For  slightly 
larger  values  of  (a   -  o    )    the  corner  Z'Z    !~ee  Fig.  7)  is 
rounded,  thus  enabling  one  to  evaluate  the  effect  of  the 
deviation  of  the  segment  PQ  from  a  sharp  corner  on  the 
deflection  angle.   In  fact,  the  results  show  that  for  a 
given  b/a  and  a/d  ratios  (see  Fig.  1)  the  effect  of  the 
curvature  is  relatively  small  as  anticipated  by  Sarpkaya  [1] 
on  the  grounds  that  the  flow  of  a  real  fluid  near  the  sharp 
corners  form  nearly  stagnant  regions  comprised  of  small 
corner  vortices  and  provides  a  slightly  rounded  streamline. 

Finally,  it  should  be  noted  that  Levi-Civita' s  method 
has  been  rendered  practicable  through  the  introduction  of 
the  Fourier-series  into  the  analysis.   It  may  be  easily 
extended  to  the  analysis  of  the  cases  where  the  lip  angle  of 
the  bucket  is  other  than  tt  or  the  bucket  has  a  cusp  on  its 
axis  of  symmetry  or  to  the  cases  where  the  impinging  jet  is 
not  free  and  emerges  from  a  prescribed  nozzle.   Such  shapes 
are  known  to  be  used  in  the  design  of  spillways,  Pelton 
wheels,  two-dimensional  cascading  thrust-reversers ,  etc. 

As  to  the  Riemann-Hilbert  method,  it  may  have  already 
been  evident  that  it  is  in  essence  a  mathematically  more 
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elegant  and  sophisticated  version  of  Levi-Civita' s  method  in 
prescribing  certain  functions.   It  is  for  this  reason  that 
no  particular  solutions  have  been  presented  and  only  a  linear 
variation  of  f(t)  was  studied.   Needless  to  say,  a  quadratic 
or  cubic  variation  of  f(t)  could  have  been  employed. 
Evidently  the  resulting  integrals  would  not  have  been 
integrable  in  closed  form.   Finally,  it  should  be  pointed 
out  that  the  Riemariii-iiilbert  method  is  perhaps  the  only 
method  aside  from  purely  numerical  methods  with  which  one 
can  incorporate  into  the  analysis  the  effect  of  body  forces 
such  as  gravity. 

In  closing  the  discussion  of  the  methods  dealing  with  the 
deflection  of  plane  jets  from  curved  obstacles,  one  must  once 
more  emphasize  that  all  of  the  foregoing  analyses  and  solu- 
tions deal  with  plane  Laplace-field  problems  where  the  real- 
fluid  effects  must  be  out  of  necessity  ignored.   Because  of 
the  limitations  imposed  upon  man's  understanding  of  turbulence, 
the  development  of  laminar  and/or  turbulent  boundary  layers 
under  a  nonzero  pressure  gradient  along  curved  walls,  entrain- 
ment  between  the  jet  and  the  surrounding  fluid  which  may  or 
may  not  have  the  same  physical  properties,  and  finally  the 
Coanda  effect  (the  tendency  of  jets  to  attach  to  adjacent 
surfaces  because  of  the  entrainment  deprivation) ,  the  real- 
fluid  effects  on  jet  deflection  cannot  as  yet  be  analyzed. 
It  is  because  of  these  reasons  that  the  analysis  of  the 
problem  with  the  assumption  of  an  inviscid  fluid  becomes 
more  important  for  1t  1s  only  through  the  comparison  of  such 
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solutions  with  the  experimantal  results  that  one  can  assess 
the  degree  of  significance  of  the  real  fluid  effects  on  the 
governing  parameters.   As  stated  in  the  introduction  and  as 
verified  by  numerous  experiments  (Robertson  [13]  >  p.  626), 
the  phenomenon  of  jet  deflection  is  governed  primarily  by 
pressure  and  inertia  forces  and  that  the  real-fluid  effects 
are  of  minor  significance  in  determining  the  forces  acting 
on  the  body  and  the  deflection  of  the  jet. .  These  facts  will 
be  stressed  again  in  connection  with  the  discussion  of  the 
analytical  and  experimental  results  for  axisymmetric  jet 
deflection  to  be  subsequently  presented. 
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III.   ANALYSIS  OF  AXISYMMETRIC  JET  IMPINGEMENT 

A.   INTRODUCTION 

As  stated  in  the  general  introduction,  the  two- 
dimensional  counterpart  of  the  jet-deflection  problem  has 
been  treated  by  several  investigators  through  the  use  of 
numerical  and  approximate  analytical  methods.   Additional 
techniques  based  on  the  modif ied-hodograph  method,  Riemann- 
Hilbert  integral  transforms,  and  the  Levi-Civita' s  method 
have  been  successfully  developed  and  employed  in  the 
preceding  section  in  the  analysis  of  the  jet  impingement 
upon  two-dimensional  curved  buckets. 

The  three-dimensional  counterpart  of  the  jet-deflection 
problem  has  not  yet  been  solved  in  any  generality.   Attempts 
to  formulate  an  exact  solution  have  been  mostly  unsuccessful 
even  for  the  axisymmetric  inviscid  flows  with  no  body  forces. 

The  case  of  a  circular  jet  striking  a  plate  normally  was 
analyzed  by  Schaach  [21]  using  approximate  methods  similar 
to  those  of  Trefftz  [22]  with  successive  adjustment  of  the 
free  streamlines  and  through  the  use  of  the  Fredholm  integral 
equations.   These  results  are  in  good  agreement  with  those 
obtained  experimentally  and  to  the  surprise  of  the  many 
investigators,  they  are  almost  identical  to  those  of  the 
two-dimensional  cases. 

Other  noteworthy  contributions  to  the  analysis  of  the 
jet-efflux  from  nozzles  and  orifices  with  straight  boundaries 
were  made  by  Southwell  and  Vaisey  [23],  Rouse  and  Abul-Fetouh 
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[24],  Garabedian  [25],  and  Hunt  [26]  through  the  use  of  the 
relaxation  and  finite  difference  techniques,  and  singular- 
integral  equations.   Schnurr  et  al .  [27]  used  the  relaxation 
method  to  analyze  the  turning  of  two-dimensional  and  axisym- 
metric  jets  from  curved  surfaces  where  there  was  only  one 
free  stream  surface;  i.e.  the  jet  was  assumed  to  leave  the 
deflector  exactly  parallel  to  the  tangent  at  the  lip  of  the 
deflector  surface  with  b  uniform  velocity  distribution  and 
the  consequences  of  the  difference  between  the  actual  deflec- 
tion angle  and  the  said  tangent  to  the  deflector  surface  was 
taken  as  a  measure  of  the  "spillage"  and  expressed  in  terms 
of  a  "turning-effectiveness  coefficient"  determined  experi- 
mentally.  Thus,  this  analysis  does  not  constitute  a  solution 
to  the  problem  under  consideration. 

Jeppson  [28]  proposed  an  inverse  formulation  of  the 
problem  for  the  flow  exiting  from  a  circular  orifice  by  using 
x  and  r  as  dependent  variables  and  cj>  and  \p   as  the  independent 
variables.   Through  the  use  of  this  technique  Jeppson  [29] 
also  analyzed  the  flow  of  a  jet  from  a  convergent  nozzle  with 
curved  boundaries. 

All  of  the  investigations  described  so  far  have  attempted 
to  provide  better  solutions  to  the  various  configurations  of 
the  orifice  flow.   These  techniques  suffer  in  general  from 
convergence  and  accuracy  problems  and  nearly  all  resort  to 
simple  trial-and-error  procedures  to  locate  the  free  surfaces 
and  to  satisfy  the  boundary  condition  that  the  free  surfaces 
be  streamlines  of  constant  velocity.   Suffice  it  to  say  that 
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a  method  was  needed  which  could  yield  solutions  of  a 
prescribed  accuracy  for  a  wide  variety  of  fluid  flow  problems 
involving  Dirichlet,  Neumann,  and  mixed  boundary  conditions. 

Zienkiewicz  and  Cheung  [30]  proposed  in  1965  the  applica- 
tion of  the  finite  element  method  to  the  solution  of  field 
problems  involving  the  equations  of  Laplace  and  Poisson. 
Since  then  a  significant  number  of  applications  of  the  finite 
element  method  to  fluid  dynamics  has  aooearsd  in  the  litera- 
ture (see  e.g.  Norrie  and  de  Vries  [31]  and  the  references 
cited  therein) . 

The  method  has  been  recently  applied  to  several  jet-efflux 
problems  involving  only  one  freestream  surface  and  relatively 
small  jet  contraction  by  Chan  and  Larock  [32].   Suffice  it  to 
say  that  the  finite  element  method  has  proven  its  versatility 
and  applicability  to  a  wide  variety  of  two-dimensional  and 
axisymmetric  flow-problems.   The  truly  three-dimensional 
problem  such  as  the  oblique  impact  of  a  circular  jet  on  a 
plane  surface,  with  or  without  gravitational  forces,  just  to 
name  one,  offers  considerable  challenges  and  demands  additional 
sophistication. 

The  special  case  of  the  gravitational  force  acting  along 
the  axis  of  symmetry  may  easily  be  analyzed  by  modifying  the 
boundary  conditions  such  that  the  velocity  along  the  free 
surface  varies  according  to  the  Bernoulli  theorem. 

This  part  of  the  investigation  encouraged  by  the  advent 
and  the  continued  success  of  the  finite  element  method  is 
devoted  to  a  determination  of  the  angle  of  deflection,  the 
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location  of  the  freestream  surfaces  and  the  velocity  and 
pressure  distributions  caused  by  a  finite  curved  boundary 
in  general  and  a  hemispherical  boundary  in  particular  placed 
symmetrically  with  respect  to  the  axis  of  an  axisymmetric , 
inviscid  jet  issuing  from  a  nozzle.   This  problem  has  not 
been  previously  analyzed  through  either  exact  or  approximate 
methods . 

As  stated  earlier,  the  deflection  of   a  jet  by  a  solid 
boundary  is  well  suited  to  potential-flow  analysis  because 
of  the  dominance  of  inertia  and  pressure  forces  in  the 
establishment  of  the  flow  pattern.   The  result  in  general 
represents  asymptotic  values  which  are  approached  as  the 
effect  of  secondary  variables  such  as  entrainment,  boundary 
layer,  compressibility,  jet  attachment  to  adjacent  surfaces 
(Coanda  effect),  etc.  are  decreased.   Specifically,  if  the 
total  angle  through  which  the  jet  is  deflected  is  determined 
for  idealized  conditions,  then  the  principle  of  impulse  and 
momentum  can  be  used  to  compute  pressures,  forces,  or  other 
dynamic  characteristics  of  the  system  such  as  the  reverse- 
thrust  ratio.   In  this  manner  one  can  not  only  optimize  the 
geometry  of  the  system  for  optimum  performance,  measured  in 
terms  of  reversed  thrust,  but  also  can  assess  the  significance 
of  the  real-fluid  effects,  cited  above,  on  its  performance. 

B.   ANALYSIS 

The  problem  may  be  defined  in  general  with  reference  to 
Figure  10.   An  axisymmetric  jet  issues  from  a  nozzle  BC 
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Figure  10   General  case  of  axi symmetric  bucket  an  nozzle 


with  an  entrance  velocity  V  at  AB .   The  jet  impinges  symmet- 


rically upon  the  curved  bucket  LKH  and  then  the  deflected  jet 
leaves  at  an  angle  $.   The  upper  and  lower  surfaces  of  the 


jet;  namely,  CDE  and  HG,  asymptotically  converge   to  the 
centerline  of  the  jet  in  accordance  with  the  conservation  of 
mass . 

It  is  apparent  from  Fig.  10  that  for  a  given  nozzle  and 
bucket  shape,  the  deflection  angle  varies  with  the  diameters 
of  the  bucket  and  the  nozzle  and  the  separation  distance 
between  the  lips  of  the  nozzle  and  the  bucket.   In  other 
words,  3  can  be  expressed  as 
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3  =  3(RQ/ro  ,  e/r  )  (56) 


Thus,  the  determination  of  this  functional  relationship 
constitutes  the  basis  of  the  problem. 

The  fluid  is  assumed  to  be  inviscid,  incompressible,  and 
free  from  body  forces.   Consequently,  the  problem  may  be 
formulated  in  terms  of  either  the  velocity  potential  <t>  or 
the  Stokes  stream  function  ifj.   Here  the  potential  function 
is  chosen  as  the  unknown  field  parameter. 

The  governing  equation  is  then  given  by 


iL!  +  Ii!+iL^=o  (57) 

a  2    r  3r    .  2    U  °'; 


in  which  x  and  r  are  the  axial  and  radial  coordinates. 
Because  of  the  axial  symmetry  there  is  no  variation  in  the 
meridional  direction.   This  equation  is  to  be  solved  with 
the  boundary  conditions  which  simply  states  that  the  normal 
component  of  velocity  be  zero  along  both  the  solid  and  free 
surfaces,  i.e.  the  free  surfaces  be  streamlines  of  constant 
velocity. 

For  an  axisymmetric  flow,  the  solution  to  the  Laplace- 
field  equation  satisfying  the  specified  normal-velocity 
boundary  condition  (8cf>/8n)   is  given  by  the  admissible 
function  <J>  which  minimizes  the  functional  (see  e.g.  Zienkiewicz 
and  Cheung  [30]) 

I(<{>)  -  Ptt   //  [(|£)2  +  (|A)2]rdxdr  -  2Ptt  6    <K|£)a  rdS     (58) 
«    dx      or  p   an 
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in  which  A  is  the  half  of  a  meridional  section  of  the  flow 
and  C  is  a  portion  of  the  curve  bounding  this  area  where  the 
normal  derivatives  are  prescribed.   The  first  integral  in 
Eq .  (58)  represents  the  kinetic  energy  of  the  fluid  within 
the  entire  control  volume  and  the  second  integral  represents 
twice  the  work  done  by  the  impulsive  pressure  p<f)  on  the 
boundaries  in  starting  the  motion  from  rest  (Milne-Thompson  [*!]) 

It  is  apparent  from  the  foregoing  that  in  the  finite- 
element  method  the  boundary  conditions  are  an  integral  part 
of  the  functional  I(<f>)  and  will  not  have  to  be  dealt  with 
separately  as  in  the  case  of  finite-difference  or  relaxation 
techniques.   This  is,  in  fact,  where  the  beauty  and  the  power 
of  the  finite  element  technique  lie . 

The  approximate  minimization  of  the  functional  I(<|>)  is 
accomplished  as  follows.   Firstly,  the  field  of  interest  is 
divided  into  M  triangular  elements  and  a  functional  I(<}>)  is 
written  as  a  sum  of  M  element-functionals  I  (<{>)  each  of  which 
has  the  form  of  Eq.  (58),  i.e. 

N 

IU)  =  I  IeU)  (59) 

1 

Figure  11  shows  the  numbering  system  of  the  corners  and  mid- 

points  and  the  corresponding  area  coordinates  £.  =  A. /A 

in  which  A.  is  the  area  of  a  subtriangle  and  Ae  is  the  area 

of  the  entire  triangle.   Then  it  is  clear  that   £,  +  Cp  +  Co  =  ^ ' 


6M 


(1/2,0,1/2) 


(1/2,    1/2,0) 


Figure    II      A  triangular  element  with   area  coordinates 


Secondly,    a   polynomial   of   sufficient    order  is    selected   to 
represent   the   variation  of  the   potential   function   (j>     within 
each  element. 

Q 

Apparently,  a  linear  variation  of  <j)   leads  to  constant 
velocities  along  a  given  line  within  an  element  and  to 
unacceptably  large  errors  in  the  velocities  between  neighboring 
elements.   A  second-order  polynomial,  however,  allows  at  least 
a  linear  variation  of  velocity  components. 

The  use  of  higher  order  polynomials  may  increase  the 
accuracy  of  the  results  and  the  convergence  of  the  minimization, 
but  it  may  also  lead  to  extremely  tedious  arithmetical  manipu- 
lations.  For  these  reasons  and  following  Felippa  [33]  and 
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Chan  and  Larock  [32],  a  second-order  polynomial  has  been 
chosen.   In  terms  of  the  area  coordinates,  this  leads  to 

6e  =  A®  X±       (i  =  1 6)        (60) 

where  A.  are  the  coefficients  of  the  polynomial  and  corresponds 

Q 

to  the  values  of  A   at  the  node  i.   It  is  important  to  note 
that  A.  must  be  regarded  as  constants  during  the  integration 
of  Eq.  (58)  in  each  element. 

The  term  x-  maY  be  written  after  Felippa  [33]  as 

xi  =  tc1(2c1-l)3C2(2c2-l)JC3(2c3"l)^C1C2^C2C3^C1c3]    (61) 

The  velocity  components  u  and  v  are  given  by 
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U  =  -r^- 

9x 

and 

v  .  M. 

dr 
Then  the  velocity  in  each  element  can  be  written  as 


ax, 

ue  =  A®  --i     (I  =  1  to  6)  (62) 

1  0  x 


3y. 
ve  -  jG 


Ae  ^-i-     (1  =  1  to  6)  (63) 

Yi  9r 


The  evaluation  of  Eqs .  (62)  and  (63)  require  the  derivatives 
of  x-  with  respect  to  x  and  r  or  the  evaluation  of  terms 
such  as 
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9x 


(H1-l)  — 1  =  (HfD  3"  =  Tl 


8x 


2A 


arU(2C0-l)  -,  3C0  b0 

*l-£ — § — 1  =  (4cP-D  -^r  ■  (4cP-i)  -^  =  t* 

8x  d  9x       d     2Ae     2 


3i_3 5 1  =  (4c  -1)  — ^  =  (4^-1)  —2-  =  T^ 


3x 


'3 


8x 


3     2Ae    ^ 


(64) 


9x 


3cP     ac 

=  I(c   — £  +  4r 


1  4 


(C-,b0  +  Cob  J  =  T 


3x    2Ae 


3(H2C3)        8C3        3C2 
8x     =  ^2  ~8x"  +  ^C3  ~9x" 


2A 


e  (C2b3  +  ^3b2}  =  T' 


8(4^^) 
3x 


Ho        3C-, 

4£ +  4c   — 

^1  8x     ^3  3x 


2A 


-  (cxb3  +  c3b1) 


in  which 


8C,    b 


l     l 


ax   2Ae 


and 


ri  -  rj 


The  foregoing  closely  follows  the  notations  developed  by 
Felippa  [33]. 


67 


It  is  evident  that  the  derivatives  of  x  with  respect  to 
r  may  easily  be  evaluated  by  changing  x  to  r  or  b.  to  a. 
given  by 


a,  =  x.  -  x. 

k    j     1 

and  denoting  the  resulting  array  of  T.  by  T. .   Thus,  one 
can  write 


(*JC   -  l)b. 

T^  =  1 1      (no  sum  on  jl)  (65a) 

1       2Ae 


Mb. c-  +  b  ,5.  ) 

T?.o  =  X   J  p  J  X  (65b) 

1  J         2Ae 


and  likewise 


4(?  -  l)a. 

T?  = (no  sum  on  i)  (66a) 

2Ae 


Tf+o  =  1  J   e  J  ±  (66b) 

itj  2Ae 


In  which 


ai      3ci 


and 


2Ae     3r 


r  =  rlCl  +  r2c2  +  1^3 
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The  minimization  of  the  functional  I(<{>)  in  accordance  with 
the  Ritz  method  requires  that  the  partial  derivatives  of 
I((J))  with  respect  to  <J>.  be  zero.   In  other  words 

M^4i=0  (67) 

which  may   be   evaluated  in  a  manner  quite   similar  to   that 
used  in  Eq.    (6*1)  .      Thus,    Eq.     (67)    reduces    to 


31    ((M    =    2P7T    //    (T?Te  +T^T!)^rdA   -    2ttP    j6    X.(|^)ardS    =    0       (68) 
S^e  ^e        1    j         1   j      j  c      1    dn 


or  to 


i^ill  =   Ke.c})e    -   Be    =    0  (69) 

9(Je  ij^j  1 


in  which  K . .  and  B.  represent,  respectively,  in  the  structural- 
mechanics  terminology,  the  element  stiffness  matrix  and  the 
load  matrix  for  a  triangular  element.   In  fluid  mechanics 
terminology,  they  simply  represent  the  geometrical  configura- 
tion of  the  grid  within  the  domain  and  the  boundary  conditions 
imposed  on  it.   These  coefficients  have  been  previously 
obtained  by  Chan  and  Larock  [32]  and  are  tabulated  in  Appendix 
G. 

In  fluid-flow  problems  where  one  or  more  parts  of  the 
boundary  are  to  be  determined  as  part  of  the  numerical 
analysis,  the  solution  of  the  system  of  equations  represented 
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by  Eq.  (69)  through  the  use  of  the  Gaussian  elimination 
procedure  requires  the  use  of  a  suitable  (i.e.  accurate  and 
convergent)  iteration  scheme.   This  important  phase  of  the 
analysis  will  be  described  in  the  following. 

C.   ITERATION  SCHEME 

The  iteration  scheme  consists  of  the  designation  of  a 
moving  and  fixed  grid,  the  adoption  of  a  scheme  for  assuming 
the  initial  position  of  the  deflected  jet,  and  the  successive 
adjustment  of  the  free  surfaces  until  the  boundary  conditions 
are  satisfied. 

1 .   Designation  of  a  Grid 

Evidently,  the  minimization  process  cannot  be  carried 
out  with  a  fixed  grid  as  it  would  be  in  a  structural  system 
subjected  to  small  deformations,  primarily  because  the 
initial  position  of  the  free  surface  is  not  a  priori  known, 
and  the  initially-assumed  boundaries  may  be  too  far  off  from 
their  final  positions.   In  view  of  the  foregoing,  the  fluid 
domain  was  divided  into  a  moving  and  a  fixed  grid  (see  Fig. 
12).   The  fixed  grid  occupies  that  region  of  the  flow  which 
is  not,  in  the  final  analysis,  to  be  intersected  by  the  free 
surfaces.   This  does  not  impose  any  restrictions  on  the  anal- 
ysis and  requires  no  more  than  a  passing  familiarity  with 
the  elementary  principles  of  jet  deflection  from  solid 
boundaries.   Should  the  boundary  of  the  fixed  grid  be  inter- 
sected by  the  free  surface  during  the  iteration  process,  a 
new  fixed  grid  has  to  be  assigned. 
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The  moving  grid  is  comprised  of  triangular  elements 
whose  nodal  points  may  be  moved  along  prescribed  lines  prior 
to  the  commencement  of  the  next  iteration.   This  procedure  may 
best  be  described  with  reference  to  Fig.  13.   At  the  end  of 
a  given  iteration,  the  nodal  points  A,  B,  C  are  moved  inward 
or  outward,  in  a  manner  and  magnitude  to  be  described  later 
and  the  positions  of  the  nodal  points  F,  G,  E,  H,  D,  and  K 
are  recalculated  such  that 

AF  =  0.15  AN 

and 

AG  =  0.50  AN 

and  similarly  for  BE,  BH,  etc.   Evidently,  the  coefficients 
0.15  and  0.50  are  arbitrary  and  different  values  could  have 
been  used  provided  that  the  resulting  elements  maintain  a 
shape  more  or  less  compatible  to  their  original  shape  and 
that  the  ratio  of  the  sides  of  a  given  element  does  not 
significantly  differ  from  unity. 

The  points  intermediate  to  AB  and  BC,  i.e.  the  points 
1  through  6,  are  calculated  at  each  iteration  through  the  use 
of  a  subroutine  called  CURVE-FIT.   This  subroutine  is  a  third 
order  polynomial  connecting  the  points  A,  B  and  C,  and  evenly 
spacing  the  points  1  through  3  between  AB  and  the  points  4 
through  6  between  BC .   It  should  be  immediately  pointed  out 
that  there  is  no  need  for  an  exotic  subroutine  to  perform 
this  procedure  and  that  any  other  judiciously-selected 
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Free  surface  and 
moving  grid  — — 


(a)   Moving  grid  along  the 
lower  free  surface 


Edge  of  the 
fixed  grid 


moving  grid  in  the 
def lected  jet 


(b)  Moving  grid  in  the  deflected  jet 


Figure  13  Construction  of  the  moving  grid 
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polynomial  could  have  been  employed.   The  sole  purpose  of 
the  foregoing  procedure  is  to  smooth  the  boundary  between 
two  successive  iterations  and  to  obtain  a  better  representa- 
tion of  the  curved  streamlines.   Furthermore,  this  procedure 
eliminates  the  occurrence  of  localized  low  or  high  velocities 
which  could  have  stemmed  from  a  zig-zag  positioning  of  the 
points  such  as  1,  2,  3  and  thereby  accelerates  the  convergence 
process.   In  the  deflected  portion  of  the  free  jet  (see  Fig. 
13b)  the  points  such  as  P  and  R  are  moved  inward  or  outward 
as  dictated  by  the  iteration  conditions  and  the  point  Q  is 
located  in  the  middle  of  PR  prior  to  each  new  iteration. 
2 .   Assumption  of  Initial  Boundaries 

The  assumption  of  an  initial  boundary  requires  the 
consideration  of  the  equation  of  continuity,  some  familiarity 
with  jet  deflection  problems  in  fluid  mechanics,  and  possibly 
the  use  of  a  deflection-angle  predictor,  if  one  is  available. 

Ordinarily,  a  deflection  angle  3  and  the  length  of 
the  deflected  free  jet  (HG)  are  assumed  (see  Fig.  14).   Then, 
through  the  use  of  V  ,  r  ,  r„,  and  the  equation  of  continuity, 
the  thickness  d  of  the  jet  is  calcuated,  to  a  first  order  of 
approximation  by 

Vr  2 

a  =  -g-2-  (70) 

Subsequently,  the  points  G  and  E  are  located  such  that  F  is 
intermediate  to  both.   Then  the  point  G  is  connected  to  H 
(the  upper  lip  of  the  jet)  by  a  straight  line,  and  the  point 
E  is  smoothly  joined  to  the  lower  portion  of  the  assumed 


7^ 


Figure  14  Nomenclature  for  a  hemispherical  bucket  and 
straight-nozzle  combination 


boundary.   In  passing,  it  should  be  noted  that  the  approximate 
use  of  the  equation  of  continuity  in  sections  intermediate  to 
GH  and  ED  may  significantly  simplify  the  assignment  of  the 
initial  position  of  ED.   For  this  purpose,  all  one  has  to  do 
is  to  assume  a  uniform  velocity  at  each  section  and  calculate 
the  thickness  of  the  jet  at  the  section  under  consideration 
through  the  use  of  an  expression  similar  to  that  given  by 
Eq.  (70). 

The  description  of  the  assignment  of  the  fluid  domain 
will  not  be  complete  without  a  brief  discussion  of  the  length 
of  the  nozzle  and  the  deflected  jet.   Clearly,  the  sections 
AB  and  EG  must  be  so  located  that  one  can  for  all  intents  and 
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purposes  assume  the  velocities  at  the  said  sections  to  be 
uniform.   Otherwise,  the  analysis  will  not  be  correct  and  the 
jet  boundary  will  be  influenced  by  the  conditions  created 
by  the  proximity  of  either  the  section  AB  and/or  EG. to  the 
points  C  and  H.   To  this  end,  preliminary  calculations  were 
carried  out  with  nozzle  lengths  (BC)  as  large  as  6r  and  it 
was  found  that  it  could  easily  be  reduced  to  3r  without 
materially  affecting  the  accuracy  of  the  velocity  at  AB .   In 
fact,  the  calculated  velocities  at  AB  with  a  nozzle  length 
of  3r  were  accurate  within  ±  0.5$.   Calculations  were  per- 
formed with  a  nozzle  length  of  3-8r  , i.e.  about  twice  the 
bucket  diameter. 

The  deflected- jet  length  GH  was  varied  between  0.5r 
and  0 . 7r   depending  more  or  less  on  the  particular  angle  of 
deflection.   Suffice  it  to  say  that  no  noticeable  change  was 
observed  in  the  calculated  position  of  the  jet.   This  is 
primarily  because  of  the  fact  that  for  a  jet  length  of  0.6r  , 
the  ratio  of  the  jet  length  to  final  jet  thickness  is  about 
3.0  and  that  according  to  all  the  previous  calculations  and 
exact  solutions  performed  for  two-  and  three-dimensional  jets, 
the  jet  attains  its  asymptotic  uniform  velocity  profile  in 
a  distance  less  than  3d. 

3 .   Iteration  of  the  Boundaries 

In  the  foregoing,  all  of  the  necessary  elements  of 
the  analysis  and  the  domain  have  been  provided  for  the  com- 
mencement of  the  iteration  process.   We  will  now  describe 
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the  process  which  eventually  leads  to  the  final  position 
of  the  jet  for  a  given  geometrical  positioning  of  a 
hemispherical  bucket  and  a  straight  nozzle  (see  Fig.  14) . 

To  begin  with,  the  free  surfaces  were  assumed  to 
be  solid  boundaries.   Thus,  the  two  boundary  conditions; 
namely,  constant  pressure  and  zero  normal  velocity  along  the 
free  surfaces,  reduces  to  a  single  condition.   Specifically, 
the  freestream  surf b.ces  -^should  be  streamlines  of  constant 
velocity.   The  normal  velocity  was  written  zero  everywhere 
along  all  the  rigid  and  free  surfaces  except  along  EG  (see 
Pig.  1*1).   An  arbitrary  value  of  4>  =  100  was  assigned  along 
EG.   Needless  to  say,  the  final  solution  does  not  depend  on 
this  particular  value  of  $ . 

To  start  a  computation  the  velocity  V   along  AB  was 
assumed  to  be  equal  to  unity  and  the  velocity  V  was  calculated 
For  an  exact  solution  V   should  be  equal  to  unity  assuming 
that  the  velocities  along  all  the  free  surfaces  are  taken 
equal  to  unity.   Evidently,  V  will  turn  out  to  be  larger 
than  unity  because  of  the  fact  that  the  pressure  along  AB  is 
larger  than  zero  and  thus  the  fluid  accelerates  along  BC  and 
decelerates  along  AL.   The  velocity  at  C  was  immediately 
corrected  to  a  first  order  of  approximation  by  iterating  on 
V  through  the  use  of  the  following  FORTRAN  STATEMENT: 

V0  =  V0  -  (VC  -  1.0)  *  V0  (7D 
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This  process  of  iteration  has  been  used  in  all  subsequent 
iterations.   Apparently,  the  iteration  on  V   could  have  been 
accelerated  by  beginning  the  analysis  with  a  V  value  smaller 
than  unity  (say,  V   =  0.9). 

Subsequently,  the  velocities  along  the  free  boundaries 
were  calculated  everywhere  through  the  use  of  the  $  values 
at  the  nodal  points  on  the  free  surface  and,  as  anticipated, 
it  was  found  that  they  were  not  everywhere  equal  to  unity. 
The  reduction  of  these  velocities  to  unity  within  a  prescribed 
error  constitutes  the  essence  of  the  iteration  process  and 
the  use  of  the  moving  grid  system. 

Let  V  represent  the  velocity  at  a  point  along  the 
free  surfaces.   Then,  the  free  surface  at  that  point  was 
moved  inward  or  outward,  depending  on  whether  V  was  greater 
or  smaller  than  unity,  by  an  amount  given  by  the  following 
statement : 


DELV  -  (VS  «*  3  -  1.0)  *  FAC  (72) 

in  which  FAC  is  a  multiplier.   The  correction  DELV  was  applied 
to  the  nodal  points  along  a  line  normal  or  nearly  normal  to 
the  boundary.   Various  values  of  FAC  have  been  used  in  the 
preliminary  calculations,  and  it  was  found  that  a  FAC  value 
smaller  than  0.05  will  yield  a  convergent  iteration.   In  the 
present  analysis  either  one  of  the  following  two  techniques 
have  been  used  in  correcting  the  free  surface.   In  the  first 
method  a  single  FAC  multiplier  of  0.015  was  applied  to  both 
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the  upper  and  the  lower  boundary,  five  iterations  were  per- 
formed, and  then  the  angle  of  deflection  $  was  increased  or 
decreased  by  two  degrees  depending  on  whether  the  velocities 
on  the  upper  free  surface  was  smaller  or  larger  than  unity. 
This  process  has  been  continued  until  everywhere  the  veloci- 
ties reached  a  value  of  1  ±  0.01.   It  is  important  to  note 
that  the  angle  of  deflection,  the  jet  thickness,  and  the 
radial  coordinate  rF  brcI  rr  change  during  each  iteration  pri- 
marily because  of  the  corrections  made  on  the  free  surfaces. 
However,  the  change  in  0    per  iteration  was  about  0.25  degrees. 
Consequently,  one  would  have  needed  at  least  40  iterations  to 
reach  a  correct  jet  position  from  an  initially  assumed  posi- 
tion with  an  error  of  ±  10°.   Evidently,  then,  the  ±  2° 
increments  applied  to  6  have  accelerated  the  convergence 
process . 

In  the  second  method  two  FAC  values  were  used;  namely, 
FACU  =  0.001  and  FACL  =  0.015-   Five  iterations  were  carried 
out  by  applying  FACU  to  the  upper  portion  of  the  jet  (along 
GH  in  Fig.  1*1)  and  FACL  was  applied  to  the  lower  boundary. 
At  the  end  of  the  fifth  iteration  the  role  of  the  FAC  values 
were  interchanged  and  five  more  iterations  were  carried  out. 
This  simply  amounted  to  first  iterating  on  the  lower  boundary 
and  then  on  the  upper  boundary.   As  before,  the  deflection 
angle  was  incremented  by  two  degrees  until  the  correct  jet 
deflection  angle  has  been  arrived  at,  as  evidenced  by  the 
uniformity  of  the  velocities  along  the  upper  and  lower  free 
surfaces. 
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The  full  description  of  the  iteration  process  would 
not  be  complete  without  stressing  few  other  computational 
procedures  which  have  enabled  the  prediction  of  the  initial 
deflection  angle  and  signaled  at  the  arrival  of  the  correct 
jet  position.   These  will  now  be  described. 

The  boundaries  were  iterated  upon  through  the  use 
of  either  one  of  the  two  methods  cited  above  until  a  deflec- 
tion angle  has  been  arrived  at  which  was  not  more  than  a  few 
degrees  off  from  the  final  jet  position.   At  that  time  no 
limitations  were  imposed  on  the  number  of  iterations  and  the 
jet  was  allowed  to  arrive  at  its  final  position.   In  passing, 
it  should  be  noted  that  the  programs  were  written  in  FORTRAN 
IV  and  double  precision  arithmetic  was  used.   Since  the  pro- 
grams on  the  Naval  Postgraduate  School  computer  (IBM  36O-67) 
had  to  be  run  overnight  for  a  period  of  about  one  hour,  the 
stops  between  two  successive  runs  have  provided  the  writer 
with  an  opportunity  to  examine  the  variation  of  the  veloci- 
ties, to  observe  the  speed  of  convergence  and  finally  to 
decide  as  to  whether  one  has  arrived  at  the  stage  of  letting 
the  computer  finally  carry  out  twenty  to  twenty-five  itera- 
tions to  establish  the  final  jet  position. 

Ordinarily,  in  dealing  with  the  analysis  of  jet  de- 
flection problems  from  arbitrarily-shaped  curved  surfaces, 
one  will  have  at  first  very  little  insight  as  to  the  value 
of  the  deflection  angle  which  the  jet  will  finally  assume. 
However,  after  one  or  two  calculations  with  one  or  two  geo- 
metrical combinations  of  a  given  bucket  and  nozzle  shape, 


80 


one  will  not  only  be  able  to  guess  fairly  accurately  the 
initial  jet-deflection  angle  but  will  also  be  able  to  derive 
a  predictor  equation.   We  will  now  describe  the  derivation 
of  such  an  approximate  equation  which  has  been  useful  in 
predicting  an  initial  deflection  angle  for  all  the  remaining 
cases  following  the  study  of  the  first  case. 

Consider  the  equation  of  momentum  applied  to  the 
fluid  domain  under  consideration: 


P   nr   2   +   iTr   2PV  2   +   frr   2pV  V.cos    g   =   T  (73) 

oo  oHo  oKoj 


in  which  P   represents  the  uniform  pressure  along  AB,  V.  the 

velocity  along  the  free  streamlines  (V.  =  1  in  the  present 

analysis),  p  the  density  of  fluid,  and  T  the  force  required 

to  hold  the  bucket  at  rest.   Evidently,  if  V  ,  3»  and  V. 

J  ?      o'   '       j 

were  known,  then  the  force  acting  on  the  bucket  could  have 
been  calculated.   In  fact,  the  present  analysis  does  enable 
one  to  do  so  not  only  through  the  use  of  the  momentum  equa- 
tion but  also  through  the  integration  of  P  dA  (dA  is  the 
projected  elemental  area  of  the  bucket  at  the  center  of  which 
the  pressure  is  P).   Since  the  velocities  along  the  bucket 
were  known  for  the  first  case,  P  dA  was  integrated  and  the 
resulting  thrust  was  calculated.   This  procedure  has  not 
only  provided  a  means  of  checking  the  thrust  calculated  from 
the  momentum  equation,  but  also  an  equation  to  predict  the 
deflection  angle  for  other  geometrical  combinations  of  the 
nozzle  and  the  bucket.   Thus,  writing 
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/  PdA  =  2  r 


v  o         v 
1  +  {~r    +  ^  cos  3 
J      J 


(74) 


one  has 


cos  3  = 


-—■  /  PdA  -  tt[1  +  (VQ/Vj)2] 


2tt(Vo/Vj.) 


(75) 


The  equations  so  far  derived  are  exact  and  Eq.  (75)  could  be 
used  as  a  mean  to  predict  the  angle  of  deflection  for  other 
R  /r  and  s/r  values  if  one  assumes  that  /  PdA  remains 
essentially  constant.   Thus  using  the  /PdA  value  obtained 
from  the  first  case  and  a  V  value  of  approximately  0.8, 
one  arrives  at  the  initial  value  of  3-   It  should  be  stressed 
that  the  foregoing  procedure  is  designed  simply  to  help  the 
analyst  to  arrive  faster  at  the  correct  solution.   It  neither 
alters  the  final  result  nor  is  absolutely  necessary  for  the 
solution  of  the  problems  discussed  here. 

D.   HEMISPHERICAL  TARGET-TYPE  THRUST  REVERSERS 

.  The  use  of  thrust  reversers  on  modern  jet  craft  for  the 
purposes  of  decreasing  the  landing  load,  maintaining  100$ 
engine  speed  while  landing,  decreasing  the  thrust  while 
maneuvering,  reducing  the  length  of  the  landing  strip,  tire 
wear  and  tear,  etc.  has  been  one  of  the  most  important 
developments  in  aircraft  industry  in  the  last  two  decades. 
To  perform  effectively,  a  thrust  reverser  should  yield  a 
relatively  small  jet-deflection  angle,  should  have  easy 
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stowage  configuration,  should  have  relatively  small  boattail 
or  stowage  drag,  and  finally  should  not  adversely  affect  the 
engine  performance.   A  hemispherical  thrust  reverser  appears 
to  fulfill  these  requirements.   Furthermore,  there  exist 
some  experimental  data  on  this  type  of  reverser  (Steffen  et 
al.  [3^]  and  Povolny  et  al .  [35])  for  comparison  with  the 
results  of  the  present  analysis. 

So  far  as  the  finite  element  method  is  concerned,  there 
is  no  restriction  on  the  nozzle  and  the  bucket  shape  and  they 
could  be  of  any  axisymmetric,  reasonably-chosen  configurations 
In  other  words,  the  finite-element  method  enables  one  to 
optimize,  within  the  assumptions  of  the  inviscid  flow  theory, 
the  relative  position  of  a  given  bucket  and  nozzle,  the  effect 
of  some  geometrical  modifications  on  the  bucket  and/or  nozzle 
shape  (such  as  placing  a  flat  plate  within  the  bucket)  on  the 
performance  of  the  reverser.   Furthermore,  it  enables  one  to 
assess  the  degree  of  significance  of  the  real-fluid  effects 
such  as  viscous  and  turbulent  energy  dissipation,  entrainment, 
Coanda  effect,  nozzle-pressure  ratio,  etc.  on  the  actual 
performance  of  the  bucket. 

The  geometrical  configurations,  namely,  the  ratios  R  /r 
and  s/r  were  chosen  partly  on  the  basis  of  the  past  experi- 
ence with  such  reversers.   For  example,  it  has  been  experi- 
mentally demonstrated  by  Steffen  et  al.  [3*0  that  the  larger 
bucket,  the  larger  is  the  reversed  thrust  and  that  the  closer 
the  nozzle  to  the  bucket,  the  larger  is  the  back-pressuring 
effect.   In  other  words,  when  a  nozzle  is  placed  very  close 
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to  the  bucket,  the  pressure  rises  in  the  nozzle  because  of 
the  blockage  of  the  flow  by  the  bucket  and  this,  in  turn, 
tends  to  affect  the  engine  performance.   Therefore,  there  is 
an  optimum  nozzle-bucket  configuration  which  yields  optimum 
reverser  performance. 

It  is  on  the  basis  of  such  considerations  that  the  nozzle- 
bucket  geometries  listed  in  Table  I  have  been  chosen  for  the 
numerical  analysis  and  experiments. 

TABLE  I 

Geometrical  Charact.         Angle  of  Deflection  3  (deg.) 

of  the  Hemispherical  s : r— r 

, ,  .     -i   j  Experimental 

buckets  analyzed      ™  ..   ,    ,    ,  ,  ,   N 
_ Finite  element  analysis (water) 

slope  of      mean     slope  of   slope  of 
the  inner    deflection  the  outer  the  outer 
Case    *V/ro      S//p   surface       angle    surface    surface 

1  1.8      0.4      21.1      20.3       19.5      18 

2  1.8      0.8      37.0      35.0      33.0      30 

3  1.6      0.8      53.2      50.5      47.8      46 


The  coordinates  of  the  free  surfaces  for  the  three  cases 
cited  in  Table  I  are  presented  in  Appendix  H.   The  computer 
program  used  in  the  study  of  the  case  number  3  is  given  in 
Appendix  I.   Suffice  it  to  say  that  similar  programs  were  used 
for  the  analysis  of  the  other  two  cases.   The  minor  differences 
between  the  three  computer  programs  essentially  resulted  from 
the  use  of  different  moving  and  fixed  grids  to  suit  the 
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particular  geometry  under  study.   An  effort  has  been  made  to 
keep  the  number  of  elements  and  the  bandwidth  in  the  resulting 
matrix  as  low  as  possible  and  yet  maintain  a  reasonable 
degree  of  accuracy.   At  various  stages  in  the  evolution  of 
the  computer  programs  larger  or  smaller  number   of  elements 
have  been  used,  and  it  was  found  that  the  number  of  elements 
used  herein  (288)  is  quite  sufficient.   Needless  to  say,  the 
coordinates  of  the  free  surfaces  may  be  calculated  to  a 
greater  degree  of  precision  by  increasing  the  number  of  the 
elements  particularly  in  the  moving  grid.   The  additional 
corrections  to  be  applied  to  the  nodal  points  on  the  free 
surfaces  at  the  end  of  the  final  iteration  were  less  than 
10   .   It  is  estimated  that  the  calculated  deflection  angle 
will  not  vary  more  than  ±1.0°.   In  view  of  the  fact  that  the 
real-fluid  effects  such  as  entrainment,  Coanda,  etc.,  can 
decrease  the  deflection  angle  by  as  much  as  ten  degrees,  it 
was  not  deemed  necessary  to  pursue  the  calculation  of  the 
deflection  angle  to  a  greater  degree  of  precision. 

A  series  of  experiments  were  conducted  with  water  through 
the  use  of  straight,  sharp-edged  nozzles  of  appropriate 
dimensions  (see  Table  I)  with  a  hemispherical  cup  of  0.72 
inch  radius  primarily  for  the  purpose  of  determining  the 
slope  of  the  outer  surface  of  the  deflected  jet.   The  slope 
was  determined  from  photographs  similar  to  that  shown  in  Fig. 
15  and  tabulated  in  column  7  of  Table  I. 

It  is  evident  from  Table  I  that  the  experimentally  deter- 
mined values  of  3  are  somewhat  smaller  than  those  obtained 
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Figure    13      >Vater   jet   deflected   by   a   hemispherical    bucket 

(R  It     -    1.6    ,    s/r     -   0.8) 
o     o  o 
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numerically  (see  columns  6  and  7 >    Table  I).   The  reasons  for 
these  differences  are  several,  the  most  notable  ones  being 
the  energy  dissipation  due  to  viscous  action  and  the  Coanda 
effect.   The  latter  will  be  discussed  in  detail  further  later 
in  connection  with  the  experiments  conducted  with  air. 
Suffice  it  to  say  that  the  Coanda  effect  with  water  is 
relatively  small  and  the  difference  between  the  computed  and 
the  measured  angles  of  deflection  is  within  the  range  of 
experimental  errors  and  numerical  approximations. 

The  theoretical  and  experimental  results  may  also  be 
compared  in  terms  of  a  measure  of  the  reversed-thrust  because 
of  its  obvious  practical  significance. 

In  the  aircraft  industry  the  performance  of  the  reverser 
is  expressed  in  terms  of  a  reverse-thrust-ratio  nR  defined  by 


n0  = 


actual  reversed  thrust 


R   forward  thrust  of  the  nozzle  alone 


or  in  other  words  by 


tt  r  2  p  y  VQ  cos  g    Vq 

u   =  ^-2 =  ^7  cos  3 

tt  rQ  p  V .  j 


The  n-n  values  found  in  the  present  analysis,  and  those 

n 

obtained  experimentally  by  Steffen  et  al .  [31!]  are  tabulated 
below. 
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TABLE  II 


Case  nR  nR 


No. 


Present  Analysis  Experimental  [34] 


0.766  0.82 

0.725  0.79 

0.515  0.62 


Evidently  the  nD  values  obtained  experimentally  by  Steffen  et 

n 

al.  [3^]  are  somewhat  larger  than  those  obtained  in  the  present 

analysis.   As  noted  earlier,  the  observed  differences  are 

primarily  due  to  the  Coanda  effect  and  the  nozzle-pressure 

ratio.   The  Coanda  effect  decreases  the  deflection  angle  and 

thereby  increases  the  reversed  thrust.   In  the  experiments 

conducted  by  Steffen  et  al .  [3^]  the  outer  surface  of  the 

nozzle  was  streamlined  in  the  form  of  a  boattail  to  decrease 

entrainment  and  thus  to  decrease  3  or  to  increase  r}~.      In 

n 

fact,  as  noted  by  Steffen  et  al .  [3^],  the  pressure  reduc- 
tions on  the  boattail  were  large  enough  to  account  for  as 
much  as  10  to  15  degrees  reduction  in  the  deflection  angle. 
Thus  the  relatively  large  differences  between  the  nR  values 
computed  in  the  present  analysis  and  those  obtained  by  Steffen 
et  al.  [34]  are  primarily  attributable  to  the  Coanda  effect. 
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In  order  to  illustrate  the  influence  of  the  Coanda  effect, 

one  particular  experiment  was  carried  out  with  air  using  a 

1.25  inch  radius  nozzle  and  a  2  inch  radius  hemispherical 

bucket  with  a  1  inch  spacing  (i.e.  with  R  /r  =1.6  and 

0  00 

s/r  =  0.8).   The  photograph  of  the  deflected  jet  (visualized 
with  aluminum  powder)  is  shown  in  Fig.  l6a.   Figure  16b  shows 
the  deflection  of  a  water  jet  from  a  similar  nozzle-bucket 
geometry.   It  is  evident  from  the  comparison  of  these  two 
photographs  that  whereas  the  slope  of  the  outer  free  surface 
with  air  is  33  degrees,  with  water  it  is  ^6  degrees.   The 
13  degree  reduction  in  the  deflection  angle  due  to  the  Coanda 
effect  corresponds  to  an  approximately  20$  reduction  in  nn  • 

n 

This  brings  the  calculated  r)„  values  into  closer  agreement 

K 

with  those  experimentally  obtained. 

The  nozzle-pressure  ratio  (total  nozzle  pressure/ambient 
pressure)  or  the  actual  velocity  in  the  nozzle  causes  varia- 
tions in  3  primarily  because  the  entrainment  needs  of  the 
deflected  jet  and  hence  the  Coanda  effect  increases  with 
increasing  jet  velocity. 

Suffice  it  to  say  that  the  analysis  of  thrust  reversers 
is  extremely  complex  for  the  cases  where  the  jet  impinges 
asymmetrically  or  where  the  gravitational  effects  are  impor- 
tant.  Furthermore,  some  of  the  practical  problems  associated 
with  thrust  reversal,  such  as  the  reattachment  of  the  jet  to 
the  nacelle  of  the  engine,  hot  gas  reingestion,  interaction 
of  the  deflected  jet  with  the  ambient  stream  and  the  adjacent 
surfaces  (in  this  case  the  outer  surface  of  the  nozzle)  cannot 
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(a)   Air  jet 


(b)   Water  jet 


Figure  16  Air  and  water  jets  deflected  by  a  hemispherical 

bucket  (R  /r  =1  .6,  s/r  =  0.8) 
o  o  o 
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yet  be  analyzed  by  any  one  of  the  existing  numerical  methods, 
even  for  relatively  simple  two-dimensional  geometries 
composed  of  straight  boundaries.   Entrainment  depends  on  the 
velocity  distribution  and  the  turbulence  level  in  the  jet, 
shape  of  the  adjacent  surfaces,  the  physical  properties  of 
the  jet  and  of  the  fluid  medium  into  which  the  jet  discharges 
Most  of  these  parameters  are  obviously  interdependent  and  at 
present  there  does  not  seem  to  be  any  hope  that  their  effect 
on  the  deflection  of  an  axisymmetric  jet  can  be  evaluated. 

Nevertheless,  approximate  analyses  such  as  the  one 
presented  herein  help  to  isolate  the  more  promising  types  of 
thrust  reversers,  investigate  their  potentialities,  evaluate 
their  ideal  performance  characteristics,  and  to  delineate 
the  range  of  importance  of  the  geometrical  variables  involved 
Then  those  cases  that  are  shown  to  be  promising  through 
computer  experiments  can  be  tested  in  a  laboratory,  at  a 
sufficiently  reduced  expense,  to  evaluate  the  effect  of 
additional  fluid-mechanical  phenomena  on  their  actual 
performance . 
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IV.   SUMMARY  OF  RESULTS  AND  CONCLUSIONS 

A  class  of  analytical  and  numerical  solutions  to  the 
problem  of  inviscid-jet  deflection  from  plane  and  axisym- 
metric  concave  surfaces  have  been  presented.   The  methods 
of  solution  for  the  two-dimensional  cases  included  the 
modified-hodograph  method  developed  herein  as  well  as  those 
devised  by  Levi-Civita  and  Riemann  and  Hilbert. 

The  modified-hodograph  method  is  capable,  within  the 
limitations  of  the  potential-flow  theory,  of  analyzing  the 
flow  in  curved  obstacles  and  generating  a  family  of  practi- 
cally suitable  bucket  shapes.   Furthermore,  this  method  has 
the  inherent  capability  of  being  extended  to  the  generation 
of  even  a  larger  family  of  concave-obstacle  shapes  through 
the  use  of  distributed  sources  and  sinks  in  the  hodograph- 
plane . 

The  method  devised  by  Levi-Civita  and  successfully 
extended  and  employed  herein  requires  the  judicious  selection 
of  the  variation  of  two  parameters  and  the  use  of  a  suitable 
series  (such  as  the  Fourier  series)  to  enable  one  to  perform 
the  integrals  leading  to  the  coordinates  of  the  obstacle  as 
well  as  of  the  free  surfaces.   A  brief  familiarity  with  this 
method  allows  one  not  only  to  generate  a  family  of  suitable 
concave  obstacles  but  also  to  round  the  corners  of  a  bucket 
composed  of  straight  segments.   This  procedure  has  been 
illustrated  by  rounding  the  corners  of  a  U-shaped  obstacle 
previously  analyzed  by  Sarpkaya  [1]. 
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As  to  the  Riemann-Hilbert  mehtod,  it  has  been  demonstrated 
that  it  is  in  essence  a  mathematically  more  elegant  and 
sophisticated  version  of  Levi-Civita' s  method  in  prescribing 
certain  functions.   It  is  for  this  reason  that  no  particular 
solutions  have  been  presented  and  only  a  linear  variation  of 
the  parametric  function  was  discussed. 

These  methods  have  provided  indirect  solutions  for  the 
shape  of  the  curved  solid  boundary  in  terms  of  a  given  jet- 
departure  and  jet-deflection  angle.   The  direct  solution  of 
the  corresponding  problems  requires  the  use  of  either  a 
sequence  of  iterations  or  a  numerical  scheme  for  solving  the 
Laplace-field  problem.   The  latter  method  has  been  reserved 
for  the  analysis  of  the  axisymmetric  cases  because  of  their 
far  reaching  and  relative  practical  significance. 

The  axisymmetric  cases,  namely,  the  deflection  of  a  jet 
issuing  from  a  nozzle  from  hemispherical  thrust  reversers 
have  been  solved  directly  through  the  use  of  the  finite 
element  method.   In  doing  so,  several  new  procedures  have 
been  introduced.   These  include  the  use  of  a  self-adjusting 
and  fixed  grid  and  a  convergent  iteration  scheme.   The  jet- 
deflection  angle,  the  slope  of  the  inner  and  outer  surfaces 
of  the  deflected  jet  as  well  as  the  coordinates  of  the  free 
streamlines  were  determined  for  a  series  of  representative 
nozzle-hemisphere  combinations.   It  has  been  demonstrated 
that  the  method  and  the  iteration  scheme  are  sufficiently 
general  to  analyze  the  more  promising  types  of  thrust-reversers , 
evaluate  their  performance  characteristics,  and  to  delineate 
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the  range  of  importance  of  the  geometrical  variables  involved. 
Finally,  a  few  experiments  have  been  carried  out  with  air 
and  water  to  assess  the  degree  of  significance  of  the  real- 
fluid  effects  on  jet  deflection.   These  experiments  as  well 
as  those  conducted  by  others  showed  that  in  cases  where  the 
Coanda  effect  is  negligible  (water  jet  in  air),  the  analytically- 
predicted  and  experimentally-measured  jet-deflection  angles 
are  in  good  agreement  and  that  the  real-fluid  effects  are 
negligible.   This  is  evidence  of  the  fact  that  the  deflection 
of  a  free  jet  by  a  solid  boundary  is  well  suited  to  potential 
flow  analysis  particularly  for  the  cases  where  the  inertia 
forces  dominate  the  establishment  of  the  flow  pattern.   In 
cases  where  the  Coanda  effect  is  not  negligible  (air  jet  in 
air) ,  the  real-fluid  effects  are  more  pronounced  and  can  cause 
the  deflection  angle  to  be  as  much  as  20$  lower  than  those 
predicted  from  the  potential-flow  theory.   Nevertheless,  the 
potential-flow  analysis  of  the  two-dimensional  as  well  as 
axisymmetric  jet-deflection  problems  enables  one  to  determine 
the  optimum  values  of  the  geometrical  variables  involved  and 
thereby  eliminates  the  use  of  costly  trial  and  error  techniques. 
With  such  information  available,  the  real-fluid  effects  could 
be  evaluated  with  a  few  carefully  chosen  experiments  and  the 
refinements  of  design  could  be  based  upon  a  secure  knowledge 
of  the  known  principles  of  hydrodynamics. 
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f  *&  *■*£  V^  ^ 

C**  APPENDIX  A                           ** 

C  #*  ** 

C**  COORDINATES  OF  THE  BUCKET  FOR  K=.C01186  COR  P  FSPGNC I NC  'l>: 

C**  TO  A  DEFLECTION  ANGLE  BETA=50.C  CEG.      ** 


X=G.05C5 
X^C.1002 
X=0.1531 
X=C2025 
X=0.2577 
X=0.3033 
X=0.3503 
X=0.4008 
X=0.45il 
X=0.5020 
X=0.5505 
X=0.6026 
X=0.6539 
X=0.7Q61 
X=0.7513 
X=0.3046 
X=0.3525 
X  =  C9095 
X=0.9545 
X=L  .0073 
X=1.0541 
X  =  l  .1092 
X=  1.1523 
X=1.2024 
X=l  .2622 
X-1.3362 
X=1.3809 
X=L .4332 
X=1.4953 
X=l  .5741 
X=l .6783 


Y=0.7916 

Y  =  l  .1097 
Y=1.3654 
Y=1.5637 
Y=1.7553 

Y  =  l  .8970 
Y=2.0305 
Y=2. L626 
Y=2.2845 
Y=2.3996 
Y=2.5025 
Y=2.6070 
Y=2.7041 
Y=2.7979 
Y=2.8753 
Y=2.9629 
Y=3.0377 
Y=3.i229 
Y=3.1875 
Y=3.2606 
Y=3.3226 
Y=3.3931 
Y=3.4463 

Y  =  3  .5059 
Y=3.5742 
Y=3.6546 
Y=3  .7011 
Y=2.7533 
Y=3.8128 
Y=3.8823 
Y=3.9634 
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C*  APPENDIX  B  ** 

C*        MODIFIED  HODOGBAPH  METHOD  ** 

C?    COMPUTATION  OF  THE  COORDINATES  OF  THE  BUCKET  ** 

C*  FOR  A  GIVEN  DEFLECTION  ANGLE  3ETA  ** 

C*  ** 

c 

IMPLICIT  REAL*8( A-P,U-Z) 

PI=3. 1415926DG 

0FI=1.0/PI 

BETA=0.0 

BETA=20.0 


C 

c 
c 


c 
c 

c 


EVALUATE  THE  BUCKET  CCGPC1NATES  EACH  10  DEGREES  OF  BETA 

DC  1  1=1,7 

BETA=BETA+10.0 
EE=BETA*PI/130.0 

EVALUATE  THE  PARAMETER  K  (XK)  FROM  EC. (17) 

XK=DEXP(-PI*CSIN(6E)/ ( 1.0-DC0S( BE)  )  ) 
WRITE(6,10) 
URITE(6,U)  BETA,XK 
1C  FORMAT (5X,  • «,/  > 

11  FGRMATt 5X, 'DEFLEC  .  ANGL.  BETA= «  ,  F5  .2 t ' P AR AM  .K= « , F 12  .  10 
0,/) 

DEFINE  XO  FRGM  EQ.(Lr>)  AND  VARY  T  FRCM  -1  TC  C 

XC=CPI*  DLCGi (1.0+XK)**2/ (4.0*XK  )  ) 

T=  1.00C01DO 

CT=C.00C1D0 

V  =  C.O 

DC  2  J=l,  1C0O0 

T=T-DT 

Y=2.*(l ./(T  +  XKJ+1  ./(T+i ./XK)-2./<T+l.))*DT/CL0G(T)+Y 

CALCULATE  AND  WRITE  X  AND  Y  COORDINATES 

X-OPI*  DL0G((T+XK)*(T+1.0/XK)/(T+1.0)**2)-XC 
WRITE (6  ,12)  X,Y 

12  FORMAT { 5X,  •X=»  , F 8 . 3  ,3 X  ,  • Y= '  ,F8.3) 

2  CONTINUE 
1   CCNTINUE 
STCP 
END 
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APPENDIX  C 

EVALUATION  OF  THE  COEFFICIENTS  OF  THE 
FOURIER  SERIES  REPRESENTING  6(a) 


Let  o(a)  be  represented  by  the  series 


6(a)  =  Z      a   cos  no 
n=0   n 


n  =  0,1,2,3,.. .      (A-l) 


According  to  Fig.  6,  6(a)  is  defined  by 


(a) 
(a) 

(a) 

(a) 

(a) 

(a) 


-  -7T 
=  IL 

2 


l°2"°l 


-  2  - 


a2"al 


2 


=  0 


tt  -a'i  i 


la2~al   a2"al 


0  <  a  <  a. 


al  <  a  <  a2 


o2  <  o  <  ? 


j    <     0     <  7T-02 


Tr-a2  <  a  <  ir-a-^ 


iT-a-]_  <  a  <  it 


(A-2) 


Multiplying  Eq.  (A-l)  by  cos  na  and  integrating  over  a  period 
one  has , 


TT 


IT 


/  8(c)  cos  na  da  =  a   /   cos   na  da 

-TT  -TT 


(A-3) 


Now  making  use  of  the  symmetry,  one  obtains 


2    7r 
a  =  —  /   6(a)  cos  no   da 
n   tt  0 


(A-4) 
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The  above  equation  can  be  written  as  the  sum  of  six  integrals 
corresponding  to  each  range  defined  by  Eq.  (A-2).   Thus,  one 
has, 


v,  "  -  \   f       (-tt)cos  no  do  +     f      5- 
n      *   I  0  a.       2 


0      -2-1 


■Vl 


a2"al 


cos  na  da 


11 
2 


ir-a. 


Tl  TT 

+  /       (-  p-)cos  na  da  +     /         (-  ^Ocos  na  da 

7T 

°2  2 


ir-a. 
+    / 


a 


TT-a     i 


■°2-°l 


°2-°l 


7T 

cos  na  da  +     /      0  cos  nada 

TT-O., 


After  the  integration  one  obtains, 


(A-5) 


a^  =  -5 [  (cos  na0 

n  (c^-a-j^ 


Xi 


cos  no^)  -  (-1)      •   (cos  nap 


cos  na-, )   ] 


(A-6) 


for  n   =    1 ,2 , 3  > •  •  • 


or 


a     = 


n  2 ,  >. 

n    (o^-o-^) 


(cos   no?   -   cos   na, ) 


(A-7) 


for  n   =    1,3,5,7 


98 


For  n  =  0,  one  has 


a  =  -  J  (A-8) 

o      2 


\ 


as  expected  from  the  average  value  of  the  function  6(a). 
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c* 

C* 

c* 
c* 

c* 


APPENDIX  D 
LEVI-CIVITA' S  METHOD 


C* 

C* 

C* 

C*: 

C 

c 
c 


c 

c 
c 


CCMf-UTATlCN  CF  THE  COCRD1AATES  OF  TFE  BUCKET  AND 
FREE  STREAMLINES,  AND  THE  DEFLECTION  ANGLE  BETA 
FOR  GIVEN  VALUES  OF  S  I  GMA1, S I GMA2  AND  K 


C 
C 

c 


c 
c 
c 


c 
c 
c 


c 
c 
c 


■.■£%Ti^%%-Jt^vi£i  s£  ;£:#;}:  ##:£:$:#  :!<  #  :jc  .-J:  ;x  s}e  #  #  :-U  v  »x  ^s  V- ~  «  re  s;:*^^^^:  5*  li  $  #  s{;  r» >;$ $  is ?;« #  $  #  ##  A 

READ    THE    VALUES    OF    SIGMA1    (SI),    SIGPA2    (S2)     />ND    K    (XK) 

READ(5,1)     S1,S2,XK 
WRITE (6, 2  J     SlfS2tXK 

1  F0RMAT(2F10.0) 

2  FCPM*T(5X,  ,SIGMA1=i  ,  F9  .  5,  4X  ,  •  S  I  GMA2=  '  ,  F9  .  5  ,  5X  ,  «K=  •  ,  F8. 

84) 

fAIN  DO-LCOP  TO  INTEGRATE   EQS . (40  ) , (41 )  ANC  (42) 

X  =  0.0 

Y  =  C.O 

PI=3. 1415926 
PI2^PI/2.0 
S  =  PI2 

DS=S/3CC.O 
DC  3  1=1,300 
S=S-DS 

LOOP    TO    EVALUATE    THE    FUNCTION    H(S),EQ.(33) 

H  =  0.0 

DC    4    J=l,2001,2 

T  =  J 

H=H+(2.0/( S2-S1) )*(C0S(T*S2)-CCS(T*S1) )#SI N ( T*S> /T**2 

4    CCNTINUE 

DEFINE  FUNCTIONS  OF  THE  INTEGRAND 

Fi=2.0*XK*XK*SIN( S)*(1.0+SIN(SJ ) / ( 1 . 0-XK*XK*CCS ( S ) **2 ) 
SMA=PI2*(S-S1)/(S2-S11 

CA=COS(  SMA) 
SA=SIN( SMA) 

SELECT    VALUES    OF    THETA    ACCORDING    TO    EQ.C39) 

IF(S.GT  .S2)  CA=0.0 

IFCS.GT.S2)  SA=1.0 

IF(S.LT.Sl)  CA=1.0 

IF(S.LT.Sl)  SA=0.0 

INTEGRATE  BY  SIMPSGN'S  RULE 

X=X+EXP(-H)*F1*CA*DS/PI 

Y  =  Y  +  EXP  (-H  )*F  1*SA*DS/P  I 
Iv  R  I  T  E  (  6  ,  6  )  S,X,Y 

6  FORMAT  (  5X,  •SIGMA=I  ,F8.4,15X,'X=,,F7.3,3X,«Y=,,F7.3) 
2  CCNTINUE 

EVALUATION  OF  THE  ANGLE  BETA 

C  =  (1.0-SSRT(1.0-XK-,;XK)  ) 

BETA=-ATAN(2.0*C/( 1.0-C-CJ ) 

DC  7  L=l,301,2 

T  =  L 

6ETA=(2./( S2-SI) ) * ( CO S ( T*S2 ) -COS ( T*S I ) ) *C** L/T**2+BETA 

7  CONTINUE 

B ETA= 18 0.0* ( 1.0+ BETA/PI) 
WRITE(6,8)BETA 

8  FCPMAT(5X, ,BETA=« , F10.3) 
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c 
c 
c 


c 
c 
c 


EVALUATION  OF  THE  COORDENATES  OF  THE  FREE  STPEANLINES 

INC=0 
R=-O.OOCl 
CR=C/1000.0 
12   CONTINUE 
X  =  C.O 

Y  =  C.O 

DO  11  1=1,1000 
R=R+DR 

F1=(XK*XK*(R*R-1  .0/R**2)/(2.0*R  )  ) /< < XK**2/4  .0)*<R+ 1 .0/ 
8R )^^2-l  .0) 

EVALUATE  W(R)  FROM  EC. (29) 

k=^ATAN(2..0*R/<  L.O-R*R)  ) 

DC  10  J=lj 1001,2 

T  =  J 

W=W+(2.G/<S2-S1) )*(C0S(T*S2)-C0S(T*SL) )*R**J/T**2 

CCNTINUE 

X=C0S(W)*F1*DR/PI+X 

Y  =  S1N(W  )*F1*DR/PI+Y 
WRITE(o,l3)  RtXtY 

FORMAT (bX,  «R  =  '  ,  F  10  .6,  5X,  «  X=  •  t  F  9  .  3  ,  2X  ,  •  Y=  ■  ,  F9  .  3  ) 
CCNTINUE 

GO  TO  12  TO  EVALUATE  THE  UPPER  FREE  STREAMLINE 

IND=IND+1 
R-l.OOOl 

DK-(C-1 e0)/1000.0 

IF(  IND. EO. 1}  GO  TO  12 

STCP 

EhC 


1C 


13 

11 
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APPENDIX  E 
EVALUATION  OF  THE  INTEGRAL  FOR  u>(t) 


M(tj  =  V-t(t-k)   /   .(a+bn)       dn        (E_1} 

71      l    (n-t)  y^nT^kT 


Noting  that 


a     bn    a  +  bt  .  , 

r  +  t  ~   z —  +  b 

n-t   n-t   n-t 


Eq .  (E-l)  can  be  written  as, 


u(t)  -  "  V-t(t-k)  {(a+bt)  /  dn       +  b  /    d     ) 

11      l       l  (n-t)  iMUn-k)      l  /  -n(n-k)  j 


(n-t)  ^/-n(n-k)     l  /  -n(n-k) 

(E-2) 


The  second  integral  in  Eq.  (E-2)  is 

k 


r      dn         -12-k  (x?    ,\ 

/   ;         =  cos    — : —  (E-3) 

1  |/n(k-n)  k 


The  first  integral  in  Eq .  (E-2)  depends  on  the  particular  range 
of  t.   For  0  <  t  <  k  ,  one  has, 
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k  d  1  I   (t+1)    ~   t   +  \fk-T  ^t(k-t) 


/      zzrzu^z:  =  —  In 


l     (n-t)  |/n(k-n)        ^  t(k-t)  |  (t-1) 

(E-4) 


For  t    <   0   and   t   >   k   one   obtains, 

k  —  -   t    +   1 

/      p=rr-   =    -  -    COS  (E-5) 

1       (n-t)  i/n(k-n)  ./  t(k-t)  t-1 


(n-t)  y/n(k-n)         V~t(k-t) 


Inserting  Eqs.    (E-3) , (E-H)    and    (E-5)    in  Eq .    (E-2)    one    finally 
obtains, 


:  (t+1)  -  t  +\jk-l  \jt(k-t) 


f.s  a  +  bt  .     2  vu  xy  V      x  V  v       '    b  iL,    .v     n    -12-k 

W(t)      = In   j- yt(k-t)       COS         — r— 

7T  _   |    (t    -    1)  *  k 


for   0    <   t    <   k,    and 

2t       .        n  

,„_*         i(a  +  bt)         -Ik                      i  b  ./ITT     7\  -12-k 

Oj(t)     =     -^ —  COS -l/t(t-k)  COS        — r— 


for  t    <   0   and  t    >   k    . 
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APPENDIX  F 
EVALUATION  OF  THE  DEFLECTION  ANGLE  3  FROM  EQ.  (52) 


Letting  t -*-«>  in  Eq .  (52)  one  has 


lim  w(t) 
t-*-°° 


i  /  (a+bnj  dn 
71  l  \/n(k-n) 


(F-l) 


and  noting   that 


3   =   iT-lim   ico(t) 


(F-2) 


Eq .  (F-l)  can  be  rewritten  after  some  rearrangement  as, 


3  =  it  + 


(a  +  —) 


k 


dn 


l  |/n(k-n) 


7,/n(k-n) 


k 


(F-3) 


Evaluating  the  integrals  and  applying  the  limits  of  integration, 
one  has, 


3  =  tt  + 


] 


/  .  bk\ 

(a  +  -p-j  cos 


-1  2  -k  ,  .  (7~ 
"IT-  +  Wk-1 


(F-4) 


Finally,  replacing  the  corresponding  values  of  a  and  b  one 
obtains, 


3  =  tt  + 


1 
2(k-l) 


,3k   n\  -1  2  -  k  , 

(^ 2)  cos     k   + 


V^  J 
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APPENDIX  G 

LIST  OF  COEFFICIENTS  K? .  AND  B® 

1J       J 


Sij    "    (aiaJ    +  biV76046 

Kn   =   3Sll(3rl   +   r2   +   r3) 

Kl2   '  K21   =  "   S12(2rl  +   2r2   +   r3} 

K13   =   K31   =  -  S13(2rl   +   r2   +   2r3) 

Kh   "  K"l   =   Sll(3rl  "   2r2   "   r3>    +   sl2<l4rl   +   3r2   +   3r3) 

K15   =  K51   =   S12(3rl  "   r2   "   2lV    +   S13(3rl  "   2r2   "   r3) 

k616  "  K6l  =  Sll(3rl  "  r2  "  2r3)   +  sl3(1"rl  +   3r2  +   3r3> 

K22   =    3S22(rl   +    3r2   +    r3) 

K?,   =   K?-    =   -   S,,(r.    +    2r,   +    2rO 
23  32  23       1  2  o 

K24  =  K42  =  S12(3rl  +  l4r2  +  3r3)  +  S22("  2rl  +  3r2  "  r3} 

K25  =  K52  =  S22("  rl  +  3r2  "  2r3}  +  S23(3rl  +  ll|r2  +  3r3} 

K26  =  K62  =  S12("  rl  +  3r2  "  2r3}  +  S23(~  2rl  +  3r2  "  r3} 

K33  =  3S33(rl  +  r2  +  3r3} 

K34  =  K^3  =  si3(_  ri  "  2r2  +  3r3)  +  S23("  2rl  "  r2  +  3r3} 

K35  =  K53  =  S23(3r'l  +  3r2  +  l4r3}  +  S33("  rl  "  2r2  +  3r>3} 

K36  =  K63  =  S13(3rl  +  3r2  +  lkr3)  +   S33(~  2rl  "  r2  +  3r3} 
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K^   =    8[S11(r1   +    3r2    +   r3)    +   S12(2r]L   +    2r2   +    rO 


+   ^(^   +   r2   +   r3)] 


K^5   =   K54   =   8s13(rl   +    3r2   +   r3} 


^S12pl    +    S22r2    +    S23r3} 


K46 


=   K 


64 


8S23(3r1  +   r2   +   rO 


MSnr1   +   S12r2   +   S13r3) 


K^   =    8[S22(r1   +   r2   +    3r    )    +   S23(r1   +    2r2   +    2r    ) 


+   S33(r1   +    3r2   +   r-   ) 


K56  =  K65  =  8S12(rl  +  r2  +  3r<3)  "  4(S13rl  +  S23r2  +  S33r3) 


K66  =  8[Sll(rl  +  r2  +  3r3}  +  Sl3(2rl  +  r2  +  2r3} 


+  S  3(3^  +  r2  +  r3) 


In  the  matrix  B.   terms,  I.    =  length  of  side  i  of  a  triangular 
element.   Thus 


Bl  =  rl 


Kdn}       *2    ^8n;   *3 


1 


B2  =  ^2 


^n;   *1    ^3n;   *3 


1 


B3  =  r3 


^n;   *1    ^n;   *2 


1 
F 
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a. 


b4        ^rl        r2ndn}      l3   3 


B 


=  (r    +  r  wi£)ai£    I 

^r2        r3M8n;      *1    3 


a. 


b6         U3  lM8n;       *2    3 
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C**  APPENDIX  H  ** 

C**  COORDINATES  OF  FREE-STREAMLINES  ** 

C**  ** 


CASE   NUMBER  ONE 


BETA=20.3 

DEG 

VO= 

=  0.817 

X  =  - 

-0.400 

k=1.000 

X=- 

-0.370 

R=1.003 

x=- 

-0.325 

R=l  .009 

x=- 

-0.30L 

R=1.012 

x=- 

-0.263 

R=1.020 

x=- 

-0.230 

R=l .030 

x=- 

-0.203 

R=1.C39 

x=- 

-0.181 

R=1.042 

x=- 

*0.  138 

R=1.060 

x=- 

-0.C99 

R=1.G30 

x=- 

-0 .066 

R=l  .100 

x=- 

-0.C38 

R=1.118 

x=- 

-0.016 

R=l.  137 

x= 

0.001 

R=1.158 

x= 

0.018 

R=l.  1.76 

x= 

0.034 

R=1.192 

x= 

0.053 

R=1.224 

x= 

0.063 

R=1.258 

x= 

0.071 

R=l  .290 

x= 

0.C72 

R=1.323 

x= 

0.072 

R=l  .J40 

x= 

0.C69 

R=l  .357 

x= 

0.C65 

R=1.375 

x= 

0.059 

R=1.392 

x  = 

0.056 

R=l  .393 

x= 

0.053 

R=  1.405 

x= 

0.050 

R=1.412 

x= 

0.  C47 

R=1.419 

X- 

0.0  44 

R=1.425 

x= 

0.040 

R=1.431 

x= 

0.036 

R=1.43  8 

x= 

0.032 

R=  1.444 

x= 

0.C28 

R=1.450 

x= 

0.023 

R=l.<»56 

x= 

0.018 

R=l  .462 

x= 

0.C13 

R=1.463 

x= 

0.007 

R=1.474 

x= 

0.001 

R=1.481 

x=- 

-0.006 

R=1.488 

x=- 

-0.013 

R=1.495 

x= 

-0.C24 

R=1.505 

x=- 

-0.036 

R=  1.514 

x= 

-0.048 

R=1.524 

x=- 

-0.C62 

R= 1.533 

x=- 

-0.07a 

R=1.544 

x= 

-0.095 

R=l  .55  5 

x=- 

-0.  113 

R=1.565 

x= 

-0.  131 

R=1.576 

x= 

-0.193 

R=1.607 

x=- 

-0.252 

R=1.634 

x= 

-0.328 

R=l  .669 

x=- 

-0.3  85 

R=1.694 

x=- 

-0.457 

R=1.724 

x= 

-0.530 

R=1.753 

x=- 

-0.6C0 

R=1.781 

x= 

-0.fc58 

R=1.805 
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x=-o 

X=-0 

x=-o, 

X  =  -0 

x=-o, 

X=-0 

x=-o 
x=-o 

X=-0 
X=-0.457 
X=-0. 528 
X=-0.5S7 


014 
C30 
C45 
061 
111 
173 
252 
310 
37  3 


R=1.805 
R=i  .807 


R=L 

R=l 
R=l 
R=l 
R=l 
R=l 
R=l 


813 
815 
83  4 
855 
884 
906 
933 


R=1.959 
R=1.934 
R=2.005 


CASE  NUMBER  TWO 


BETA=25.0 

DEG 

VO=0.£85 

X=-C. TOO 

R=l  .000 

X=-0.774 

R=1.C01 

X=-0.743 

R=l  .003 

X=-0.723 

R=1.006 

X=-0.697 

R=1.C09 

X=-0.663 

R=1.013 

X=-0.629 

R=1.C19 

X=-0.595 

R=1.025 

X=-0.561 

R=1.030 

X=-0.525 

R=1.038 

X--0.489 

R=l  .047 

X=-0.453 

R=1.057 

X=-0.418 

R=1.069 

X=-0.394 

R=1.080 

X=-0.270 

R=1.C89 

X=-0.346 

R=l  .100 

X=-0.324 

R=ielll 

X=-0.301 

R=L.  125 

X=-0.281 

R=l  .139 

X=-C260 

R=1.153 

X=-0.241 

R=1.168 

X=-0.224 

R=1.102 

X=-0.208 

R=1.199 

X=-0.  194 

R=1.216 

X  =  -0..  181 

R=1.235 

X  =  -0.  166 

R-l.253 

X=-0.155 

R=l  .282 

X=-0. 146 

R=1.308 

X=-0.  142 

R=1.335 

X=-0.140 

R=1.358 

X=-0.  140 

R=1.382 

X=-0.143 

R=  1.404 

X  =  -0.  149 

R=1.425 

X=-0.  154 

R=1.443 

X=-0.161 

R=l  .461 

X=-0.  168 

R=1.478 

X  =  -0.  177 

R=1.495 

X  =  -0.  184 

R=l  .506 

X=-0.  191 

R=1.517 

X=-0.200 

R=1.529 

X=-0.209 

R=1.539 

X=-0.243 

R=1.579 

X=-0.304 

R=l  .635 

X=-0.367 

R=1.689 

X=-0.430 

R=1.733 

X=-0.494 

R=1.791 

X=-0.562 

R=1.838 

X=-O.COO 

R=l  .800 

X=-0.068 

R=1.8C0 

X=-0.083 

R=1.803 

X  =  -0.  137 

R=1.809 

X=-0.  150 

R=1.855 
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X=-0.234  R=1.883 

X=-0.288  R=1.947 

X=-0.367  R=1.976 

X=-0.425  R=2.037 


CASE  NUMEER  THREE 


6ETA=50.5 

DEG 

V0=0.810 

X=-0.800 

R=l 

..000 

X=-0.765 

R=3 

.  .007 

X=-0.729 

R=] 

..015 

X=-0.694 

R=] 

..025 

X=-0.660 

R=J 

.  .036 

X=-0.622 

R=l 

..045 

X=-0.603 

R=] 

..055 

X=-0.575 

R=] 

.  .066 

X=-0.543 

R=i 

..078 

X=-0.525 

R=  ] 

.  .090 

X=-0.5C3 

R=] 

..102 

X=-0.482 

R=] 

..115 

X=-0.46i 

R=] 

.  .130 

X=-0.451 

R=] 

..138 

X=-0.440 

R=J 

..147 

X=-0.430 

R=] 

..155 

X=-0.420 

R=j 

..165 

X=-0.412 

R=] 

1.173 

X=-0.4C4 

R=i 

..182 

X=-0.296 

R=l 

..191 

X=-0.388 

R=i 

.  .200 

X=- 0.3  84 

R=] 

..206 

X=-0.379 

R=] 

L.212 

X=-0.375 

R=l 

..219 

X=-0.37l 

R^] 

..225 

X=-0.366 

R=J 

L.234 

X=-0.361 

R=J 

..242 

X=-0.257 

D  — 

..257 

X=-0.353 

R=] 

.  .260 

X=-0.349 

R=l 

..268 

X=-0.346 

R=] 

L.27  6 

X=-0.343 

R=] 

..284 

X=-0.34  1 

R=] 

..293 

X=-0.339 

R=i 

l.300 

X=-0.337 

R=] 

.  .3U8 

X=-0.235 

R  = 

L.316 

X=-0.334 

n  — 

IN  —    J 

l  .325 

X=-0.333 

R=] 

..331 

X=-0.232 

R=. 

L.337 

X=-0.331 

R=] 

L.344 

X=-0.331 

R=J 

L  .350 

X=-0.327 

R=j 

L.443 

X=-0.367 

R=] 

L.513 

X=-0.399 

R=j 

L.cll 

X=-0.495 

R=] 

L.742 

X=-0.561 

R=] 

L.830 

X=-0.627 

R=. 

L.918 

X=-0.000 

R=j 

L.60'0 

X=-C.C90 

r=; 

L.58  8 

X=-0.101 

r=: 

L.654 

X=-0.  158 

R  = 

L.678 

X=-0.175 

R  = 

L.734 

X=-0.255 

R  = 

L.809 

X=-0. 122 

R  = 

L.697 

X=-0.397 

R  = 

L.977 

X=-0.469 

r=; 

'.059 

110 


Qif.  A  ##*  ^^^i:**^^:^*?**:*::^:}:***:*-^  Y.«#A*:fc#:?ici«:fc:*A##:=:fc>!«#:£:$:0#»*:?::it>it*:*# 


c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 


APPENDIX     I 

FINITE    ELEMENT    METHOD 

COMPUTATION  OF  THE  ANGLE  CF  DEFLECTICN  ANC  THE 
CCCPDINATES  OF  THE  FREE  STREAMLINES  FOR  AN 
HEMISPHERICAL  BUCKET 


** 
** 
** 

** 
** 
** 
*# 

#* 


c 

IMPLICIT  REAL*8< A-H,0-Z) 


C 
C 


DINENSICN  ING(288,6),JET(9,3),JABAU2,6),X(6  83),R<683) 
8,RR(65)  ,XX(65)  ,ELE(6,6),Z(683,6G)  ,B(683)  , AN  OLE (45) 
8,V(80) 

VALUES    FOR    THE    NUMBER    OF    ELEMENTS     (  NUMEL  >  ,  B  ANDMDTH     (NBANI 
NUMBER    CF    NODAL    POINTS     (NNODE),NODE    AT    LIP    OF    ELCKET     (N2) 

NKCDE=683 

NEAN=60 

N2=597 

MNEL=2S8 

P 1=3. 14L5926D0 

DK  =  1.0  DC- 
ASSIGN    STARTING    VALUES    FOR    DEFLECTION    ANGLE     (BETA)    AND 
VELOCITY    AT    THE    ENTRANCE    OF    THE    NOZZLE    (VO) 

BETA=0. 8CDC 

VC=0.85C0 

FACU=0.001D0 

FACL=0.C2DQ 

BEAD    EQUIVALENCE    CF    NODAL     POINTS 

READ(5,1)     (     I , ( INO( I, J)  ,J  =  1,6)  ,N=1,NUMEL) 

1  FCPMAT (3(713) ) 

READ    COORDINATES    OF     INTERIOR    POINTS 
READ(5,2)< I,X(I),R(I) ,N=1,90) 

2  FCRMAT(  5(  13  ,F6.0, F5  .0  J  ) 

READ    X    VALUES    OF    THE    NODES    OVER    THE    BUCKET 
READ(5,3)     (     I ,X(I  )  ,N  =  1  ,18) 

3  FCRMAT(9(  I3tF5.0)  ) 

READ    EQUIVALENCE    CF    NCDES    IN    THE    JET    SECTION 
READ(5,4)-    (<JET( I , J) , J=l,3) ,1=1,9) 

4  FCRMAT(8(3I3) ) 

READ    ECUIVALENCE    OF    NODES     IN    THEFLOATING    GRID 
READ(5,5>( (JABA(I,J),J=1,6),I=1,12) 

5  FORMAT ( 4( 613) ) 


C    READ    ThE    COORDINATES    CF    THE    ASSUMED    FREE    SURFACE 
READ(5,6)  (I,X(  I ) ,R(I)  ,N=1,10) 
6   FCRMAT (5( I3,F6.0, F5.0) ) 


LOCATE    THE    NODES    OVER    THE    EUCKET 

RADI0=1 .6000 

DC  40  1=63,78 

J=INC( I ,3) 
4  0       RC  J)=DSCRT(RADI0**2-X(J)**2  ) 

LOCATE    THE    UPPER    FREE    SURFACE    REALTIVE    TO    POINT    N2 
DX=0.G77D0 
DISP=0. 0 
DC     10     1=55,59 
XC=61-I 

K  =  INQ( I  «  1) 
X(K)=X(N2)-CX*XD 
1C   R(K)=R( N2)-DTAN(6ETA)*(X(K)-X(N2)+CISP) 
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c 

C  LCCATE  THE  PCINTS  CLOSE  TO  THE  LIP 

ZR=R(647)-R(N2) 

X(617)=X(647)*0.25D0 

X(fc27)=X(647)*O.5D0 

X(637)=X(647)*0.75D0 

P(617)=R(N2)+ZR*0.25D0 

R(627)=R(N2)+ZR*0.50D0 

R(637)=R(N2)+ZR*0.75D0 

READ(5,6)U,X(I)  ,RU)  ,N=lt16) 

RAL  =  9 
C 

C  STARTING  POINT  FOR  ITERATIONS  ON  (BETA)  AND  (VC) 
10CC  CONTINUE 

INCEX=0 
C 
C  LOCATE  THE  REGION  C  OF  THE  LOWER  FREE  SURFACE 

TPES=0.05D0 

VETA=BETA+TRES 

DC  20  1=1, MAL 

I1  =  JET(  1,1  ) 

I2=JET(  I  ,2) 

I2=JET(  1,3) 

IHI.GE.8)  GO  TO  21 

R(  I2)=DSGRT(R(I1)**2-VC*DC0S( VETA) ) 

X(I2)=X(I1)-DTAN(VETA)*(R(I1)-R(I2)) 
21   X(I3)=(X(J 1)+X( 12)  )*0.5D0 
2C   RU2)  =  (RII1)+R(  12)  >*C.5D0 

X(591 )=C.06D0 

R(5S1)=1.48D0 

X(606)=-0.  16D0 

R(606)=l .48D0 


c 
c 

c 

L 

IF(MAL.EU.S)  X(618)  =-0«.33D0 
IF(MAL.EG.9)  R(6  18  )= L.44D0 

STARTING  POINT  TO  ITERATE  THE  FREE  SURFACES 
100  CONTINUE 

INCEX=INDEX+1 

LCCATE  THE  INTERIOR  POINTS  IN  THE  FLOATING 
DC  30  1-2,12 
I  i  =  JA£JA(I  ,  1) 
I2=JABA(I ,2) 
I2=JABA(I  ,3) 
I4=JAbA(I ,4) 
I5=JABA(I,5) 
I6  =  JABA(I  ,6) 

GRID 


X(I2)=X(UH0.15D0*(X(I4)-X(I1)  ) 

X(  13)  =X  (II  )+0.5D0*(X(  I4)-X( II) ) 
R(I2)=R(I1)+0.15D0*(R(I4)-R(I1) ) 

.  R(I3)=R(I1)+0.50D0*(R(I4)-R(I1)  ) 

X(I5)=C.5DG*(X( I6)+X( 12)  ) 
3C   R(I5)  =  C.5DC;MR(I6)+R(  12)) 
C 
C  LCCATE  INTERMEDIATE  PCINTS  IN  ASSUMED  LOWER  FREE  SURFACE 

DC  50  1  =  1,  1.2 

K  =  JABA(  1,1) 

X  >  ( I )  =  X  ( K ) 
5C   RR(I)=R(K) 

CALL  DRVF IT ( 2 , 12 , RR ,XX , 4 , 45 , RR , XX, ANGLE ) 

DC  60  1=9,48 

K=  INCH  I  ,3) 

X(K)=XX (1-3) 
6C   R(K)=RR( 1-3) 

IF(KAL  .NE..9)  GO  TO  28 

X (fcl8) =-0.33D0 

R(618)=1.44D0 

X (628)=-0.36D0 

R(628) =1.52D0 

X(638)=-0.41D0 

R  (638)  =  1.60D0 
28   CCNT1NUE 
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RAO57.2958D0 

ALFA3=DATAN( ( R( 6  72 ) -R ( 647 )  ) /( X ( 647  )-X ( 672 )  )  ) *RAD 
WPITE(6,67)     ALFA3,ALFA3 

CALCULATION1    OF    THE    DEFLECTION    ANGLE     ALFA 

ALFA1=DATAN((R(672)-R(667))/(X(667)-X(672)))*RAD 
ALFA2=DATAN<  (  R  (  677  )  -R  (  657  )  )  /  (  X  (  6  57  )-X(  677  )  )  )*RAD 
GAFA1=DATAN( ( R ( 668 ) -R ( 64tf )  ) / ( X ( 6  4 8  )  -X ( 66 8  )  )  )  * R AO 
GAMA2=DATAN(  (R(658)-R(638))/(X(63£)-X(658))  )*RAD 
WfUTE(6,61)  ALFA1,GAMA1 
hRITE(6,61)    ALFA2,GAMA2 

6L       FCRMAT(5X,'    UPPER    FSL  •  ,  F  10  *  2  ,  5X  ,  «     LOWER    FSLSF10.2) 
GAPA3=DATAN((R(666)-R(658))/(X(658)-X(668)))»RAD 
GAMA4=DATAN(  (  R(  65  8  )-R  (  64  8)  )  /  (  X(  64  8  )-X  (  6  5  8  ))  )  *R  AD 
WRITE (6,67)     GAMA3,GAMA4 

67       FCRKA'f  (2P2C.2) 


KPITE(6,500)  VOtBETA 
5CC  FCRMATt /,5X,  '  V0=  •  ,F  J  0.  4  ,  5X  ,  ■  BETA  =  •  ,  Fl  0  .  4  ,/  ) 


ZERO  CUT  MATRICES  (Z) 
DC  42  L=L,NNODE 
DC  43  M-l.NBAN 
43  Z(L,M)=0.0C0 
42  B(L)   =C.CDO 


EVALUA 
DC 
I  1 
12 
12 
Al 
A  2 
A3 
E  L 
B2 
B2 
AR 
SI 
Si 
SI 
S2 
S2 
S3 
EL 
EL 
EL 
EL 
EL 
EL 

81  1 
EL 
EL 

8)- 
EL 
EL 

8(  1 
EL 
EL 
EL 
EL 
EL 

82. 
EL 
EL 

811 
EL 
EL 

8R( 
EL 
EL 
EL 


TE  T 

44 
=  1N0 
=  INO 
=  INO 
=  X(  I 
=  X(  I 
=  X(  I 
=  R(I 
=  R(  I 
=  R(  I 
EA=0 
1=(A 
2=  (A 
3=  (  A 
2=(A 
3=  (A 
3  =  (A 
E(  1 
E(l 
E(2 
E(  I 
E(3 
E(  1 
)+3 
E(4 
E(  1 
2.0 
E(5 
E(  L 
l)+3 
E(6, 
E(2, 
E(2, 
2(3, 
E(2, 
0*R( 
E(4, 
E(2, 
)  +  14 
E(5, 
E(2, 
11)  + 
E(6, 
E(3, 
E(3, 


HE  E 
K  =  l, 
(K,  I 
(K,2 
(  Kf3 
3)~X 

1  )-X 
2)-X 

2  )  -  R 
3)-R 
1  )-R 
.5~M 
1*A1 
1  ■■<  A  2 
1*A3 
2*  A  2 
2*  A3 
3*  A3 
1  J  =  3 
2)=- 
1  )  =  E 
3)=.~ 

1)  =  E 
4)=S 
0*R( 
L)  =  E 
5)=S 
R  ( 1 2 
1  )  =  E 
6)  =  S 
.C*R 
1)=E 
2)=3 
3)=- 

2)  =  E 

4)  =  S 
I  1)  + 
2)  =  E 

5)  =  S 
.  C#  R 
2)  =  E 
6)^S 
3.0:' 
21  =  ^ 
^>)  =  3 
4)  =  S 


LEMENTS 

NUMEL 

) 

) 

) 

(12) 

(13) 

(III 

(13) 

(ID 

(12) 

A2*B1-A 

+B1-BL ) 

+B1*B2) 

+Bi*B3 ) 

+B2  *B2) 

+B2*B3  ) 

+B3*B3  ) 

.0*S11* 

SL2*(2. 

LE(1,2) 

S13*(2. 

LE(1,3) 

LI* (3.0 

12) +3.0 

L  E  (  1 , 4  ) 

12i:(  3.  •-* 

)-R( 13) 

LEU, 5) 

11* (3.0 

(12) +3. 

LE(L ,6) 

.0*S22^ 

S23*(R( 

LE(2,3) 

L2*(3.0 

3.0*R(  I 

LE(2,4) 

22*(-R( 

(12) +3. 

L  E  (  2  ,  5  ) 

12* (-R( 

U  (  I  2  )  -  R 

LE(2,6) 

.>:533*M 

13*(-R< 


AND  (B) 


OF  MATRIX  (Z) 


1  *  B  2  ) 
/  (ARE 
/(ARE 
/  (ARE 
/(ARE 
/(ARE 
/  (ARE 
(  3 .  0* 
C*R(  I 


A-60.0) 

A*60.0) 

A*60.0) 

A*60.0) 

A*60.C) 

A *6 0.0) 

R(I 1J+R(I2)+R(I3) ) 

1)+2.0*R(  I2)+R( 13) ) 


0*R(  I  l)+R(  I2)+2.C*R(  13)  ) 

»R(I1)-2.*R(I2)-R(13))+S12*(14.G*R( 
*R(12) ) 

R(I1)-R(I2)-2.*R(I3))+S13*  (3.0*R(  IL 
) 

*R(  Il)-R( 12)-2.C*R(I3)  )+S13*( L4.0*R 
0*R(  13)  ) 

( R< I l)+3.0*R( I2)+R( 13) ) 
U)+2.0*R(  I2)+2.C^R(  13)  ) 

*R(I L)+14.0*R( I2)+3.0*R( 12) )+S22*(- 

2  )-R(  13) ) 

I  1)+3.0*R( I2)-2.0*R(I3)  )+S23*(3.*R( 

0  *  R  ( 1 3  )  ) 

1  1)+3.0*R(  I2)-2.C*R(I3)  )+S23*(-2.0* 
(13)  ) 

R(  Il)+R(  1'2)+3.0*R<  13)  ) 

I  1)-2.0*R(  I2)+3.0*R(I3)  )+S23*(~2.0* 


113 


8R( 
EL 
EL 

8R  ( 
EL 
EL 

82. 
EL 
EL 

8(1 
EL 

81  ) 
EL 
EL 

61) 
EL 
EL 

8+2 
EL 

611 
EL 
EL 

6(  I 


Ill-R 
E  (  4 ,  3 
E(3,D 
J.I) -2 
E(5,3 
E{3,6 
0*R(  I 
E(6t3 
E(4,  4 
l)+2. 
E(4,5 
+  S22  < 
E  (  5  ,  4 
E  ( 4  ,  6 
+  S12- 
E(6,4 
E(5,5 
.0*R( 
E  (  5 ,  6 
+S23* 
E  ( 6  v  5 
E(6,6 
1)  +  R( 


(12 
l=E 
)  =  S 

•o* 

)  =  E 
)  =  S 
1)- 

>  =  E 
)  =  8 
0*R 
1  =  8 

R(I 
)  =  E 

)=B 
P(I 
)=E 
)  =  8 
12) 
)  =  8 
R{  I 
)  =  £ 
)  =  8 
12) 


)+3.0*R( 
LE(3,4) 
23*(3.0* 
f<(  12  J  +3  . 

LT (3,5) 
L3*(3.0* 
R(  12) +3. 

LE(3,6) 
.0*(SU* 
(  I2)+R( I 
.0*S13"M 
2)+S23*R 
LE(4,5) 
.0*S23*( 
2)+S13*R 
L  E  (  4  f  6  ) 
.0*1 $22* 
+2.0*R< I 
,0*S12*< 
2)+S33*R 
L  E  (  5  ,  6  ) 
.0*( Sll« 
+2.*R(I3 


13)  ) 

R( I l)+3.0*R(I 2)+I4.0*R ( 12) )+S33*(- 

0*R<  13) ) 

R( I l)+3.0*R( I2)+14.0*R( 13 ))+S33*(- 

0*R( 13) ) 

(R(  I  1  )+3.0*R(  I2)+R( 13) )+S12*(2.0*R 
3)  )+S22^(3.*R(  Ii  )+R( I2)+R( 13) )  ) 
R(I  1)+3.0*R( I2)+R( 13)  )-4.*(S12*R(I 
(13)  ) 

3.0*R(I l)+R(I2)+R(I3) )-4  .*( Sli*R( I 

(13)  ) 

(R(I1)+R(I2)+2.C*R(I3))+S23*(R(I1) 
3)+S33*(R(  11  )+3.*R(  12  )+R(  13)  )  ) 
ft( I  1)+R( I2)+3.C*R( 13)  )-4.C*(S13*R( 
(13)) 

(R(  I  l)+R(I2)  +  2.C*R( 13)  )  +  S13*(2.0*R 
l)+S33*(3.*R(  Il)  +  R( I2)+R(  13))  ) 


C 

C 


LC£D  MATRIX  (Z) 

DC  45  J =1,6 

11  =  INCH  K,J) 

DC  46  L-1,6 

III=IN0(K,L) 

X  F I  I I.GT.I  ii )  GO  TC  46 

NW=III-II+1 

IF(NW.GT.NBAN)  N6AN=NW 

Z( II,NW)=Z(I I,NW)+£LE( J,L) 
46  CONTINUE 
hi    CONTINUE 
kk    CONTINUE 

LC*D  VECTOR  (B)  WITH  BOUNDARY  CONDITICNS 

AT  THE  INLET  OF  THE  NOZZLE 
B(l)=-0.0 

B(2)=-R(3)*R(3)*V0/3.C 
E(3)=-R(3)*R(5)*V0/6.0 

B(4)=-(  R(  S)  +  R(3)  )*{R(5)-R(3)  )>:VC/3.0 
B(5)=-rR(5)  *(R(7)-R(3)  )*V0/6.G 
B(6)=-(R(.7)+R(5)  )*<R(7)-R(5))*VO/3.0 
B(  7)=-R<7)~HR(9)~R(5)  )*VQ/6.0 
B(8)=~(R(9)+R(7) )* (R(9)-R(7) )*VO/3.0 
B(9)=-R(9)*(R(9)-R(71)*V0/6.0 

AT  TEE  EXIT  OF  THE  JET  (TAKE  POTENT  I AL=  100. 0 ) 

DC  70  1=668,672 

B( I)=10.0D12 
70   Z (1 , 1)=  10.CD11 

SOLUTICN  OF  THE  POTENTIAL  EY  GAUSSIAN  ELIMINATION 
NAL  =  2 
CALL  BANDEC(Z,B,NNCCE,NEAN, 1 ,683,60,  1) 

CALCULATIONS  FOR  THE  VELOCITIES  ALONG  THE  BOUNDARIES 

DC  110  1=1,79 

J  1=  INCH  I  ,  1) 

J3  =  IN0(  I  ,3) 

S=DSQRT<  (R(  Jl)-R(  J3)  )**2+(X(  J1J-XU3)  )**2) 

V(  I  )=(B( J3)-B( J  1)  )/S 

IF(I.GE.55)  V(I)=-V(I) 

IF(I.GF.55)  J3=J1 

WRITE(6,ill)  J3,X(J2),R(J3),V(I) 
111  FCPPAT(5X, 'NODE'  ,  15,5X,«  X=',  F 7  .3 , 3X , r R=  •  ,  F7.3,10X,« 

8  VELOC=« ,F7.3,/) 
110  CONTINUE 


ll'l 


c 
c 

c 
c 


c 
c 


c 

c 


VELOC.  AT  LIP  CF  NOZZLE  (VI)  AND  LIP  OF  BUCKET  <V2) 
V1  =  (BU94)-BQ99)  J/0.0475D0 

V2=V<63) 

WFITE(fc,U2)  VI,  V2 
112  FGKMAT(//»  '  VELOC.  L  I  P 1  =  •  , F 10 . 4 , 5X  ,  « 
8F1C.4,//) 


VELCC.  LIP2=» 


CORRECTION  FCR  THE  UPPER  AND  LOWER  FREE  SURFACES 
WRITE(6, L22) 

FCP  TFE  LOWER  FREE  SURFACE 
DC  120  1=9,54 
AN=PI-oETA 

1FU.LE.48)  AN  =  ANGLE(  1-4) 
K  =  I  NO  (  I  ,  3 ) 

DELV=(V  (I  )**3-1.0)*FACL 
WR1TE(6,121)  K,DELV 
XX  (I  )=X(K)-DELV*DSIN(AN) 

120  RR<I)=R(K)+DELV*DCOS(AN) 

FCR  TFE  UPPER  FREE   SURFACE 

DC  L30  1=55,62 

K=INO(  I  ,1) 

DELV=(V<I )**3-1.0)*FACU 

WRITE(6,12i)  K,DELV 

XX  ( I)=X<K)+DELV*DSIN(AN) 
130  RR U)=R(K)-OELV"DCCS(AN) 
122  FCRMAT(/,»  CORRECTIONS  TO 

121  F0RMAT(5X,«  NODE  •  , I  5 , 5 X  , ' 
VC=V0-(V1-1.0)*V0 


BE  MACE  AT  EACH 
CORRECTICN1 ,F8. 


NODE 
4) 


/) 


OF  THE 
1=9,62 


RELOCATION 

DC  140 
(</  =  ra 

IF{I.GE.55)  M=L 
K=INO(I fM) 

X(K)  =  XX(  I  ) 
R(K)  =  RR(  I  ) 
140   CCNTINUE 

DC  160  I=8C,99 
I  1=  INCH  1,1) 
I2-IN0(  1,2  ) 
I  3= I  NO (  I 
IFd.GE 


NEW  FREE  SURFACE 


,3) 

94) 


I 


161 

lee 


13 

V3=(6(  I2)-B(  U))/C 
WRITE (6,  16  1)  12, V3 


(  12)-X(  ID  ) 


CONTINUE 


«VELOC=l  ,F9.3) 


CORRECT  VO 
GO  TO  101 


AND  THE  ANGLE  EETA 


EACH  FIVE  ITERATION 

IF( INDEX. EC. 10) 

IFi INDEX. NE. 5)  GO  TO  100 

FACU=Q.015D0 

FACL=0.001DO 

GC  TO  100 
101  CCNTINUE 

FACU=0.C01D0 

F/>CL=0.C15D0 

CORRECTION    CF    THE    ANGLE    BETA 

DEE=0.0300 

BETA=BETA+DBE 

DC  19  1=55,59 

XC=61-I 

K=1N0(I  ,1) 

X(K)=X(N2)-DX*XD 
19   R  (rx)  =  R(N2)-L,TAN(3ETA)*(X(K)-X(N2)  +  DISP) 

ZR=R(647)-R(N2) 

X(£17)=X(647)*0.25DO 

X(627 l=X(647)*0.5D0 

X(637)=X(647)*0. 75D0 

R(617)=R(N2)+ZR*0.25D0 


.  115 


P(627)  =  MN2)  +  2R*0.S0D0 
R(637)=R(N2)+2k*0.  7  5D0 
GC    TO    LOOO 
SSSS    CONTINUE 
STOP 
Ef\C 
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10 


40 


IOC 
199 


19C 
200 

10CC 
1010 


I 8* CUT 
fFLICI 
I  MENS  I 
CCF=NE 
C  100 
r(DABS 
=  1  +  1 
=  MNO( 
C  100 
=J+2-M 
=A(I,L 
C  40  M 
(J,  MM) 
N=MINO 
C  100 
K=L+K- 
(  J  t  K  )  = 
C  199 
(KEG, I 
C  200 
=NEC-I 
=  MNO( 
C  200 
C  190 
=J+L-1 
(J,  MM) 
(J,  MM) 
ETURN 
k  R I T  E  { 
FCRMAT 
!\D 


INE  BANDECf A,F,NEG,MAXB,NVEC»NROWfNGAN,NBEC) 

T  REAL*8( A-r,0-Z  ) 

CN  A(MROWtNEAN)  ,  F(NR0W,N8EC) 

L'-l 

1=1, LOOP 

(A(  1,1) J-l.E-30  )  1000,1000,10 

I+MAXB-1.NEC1 
J  =  M,N 

)/A(Ivl) 

f'=l,NVEC 

=  F( J,MM)-D*F(  I, MM) 

(MAXB-L+1 ,NEC-J+l) 

K=1,MM 

1 

A( J,K)-D*A(  I,NN) 

I=1,NVEC 

)=F(NEC, I )/A(NEG,l) 

I=2»NEQ 

+  1 

KEC-J+1,MAXE) 

PM=1,NVEC 

L  =  2,K 

=F( J,MM)-A( J,L)*F(M,MM) 
=F(J,MM)/A(J,1J 

6,1010)  I 

(29H1   ZERO  CIAGONAL  ELEMENT  ROW     14) 
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SLEROUT  INE   DRVF  IT ( N C, L ,X , Y ,M ,N , U , V, ANGL ) 


I 
C 

8, 
E 

1 
R 
E 

1 
C  PR  EL  I 
1C  f 
M 
L 
L 
M 

N 
I 
I 
I 
I 
I 
G 
I 
I 
12  D 


NPLICIT  REAL*8( 
INENSICN     X(L 
ANGL(N) 
CUIVALENCE 


EAL*8 

GOIVALENCE 


BINARY 
CC  =  MD 
DM  =  MD 
C  =  L 

Nl=LO- 
C=M 

P i-MO- 
C  =  N 


(Ml 
(X2 

Mir 

(W2 

(A2 
PRCCESSI 


A-H,0~Z) 
),Y(L)fU(N) ,V(NJ 

,Bl),(M2,B2),(M2,E2),(M4,B4), 

,PO)  ,  (Y2  ,QO) ,(T2,G1) 

M2,N3,f4 

,Q2),(W3,Q3),(A1,P2),(B1,P3), 

,  D  Z )  ,  (  S  U  ,  R  ,  Z  ) 

NG 


I  1 


12 

G 
14  C 

15 

G 
16  D 


17 

18 


IS 
R 
R 
C  MAIN 

ac  k 

c 

C  RCLTI 
C 

2C 


21 


4C 


(MDO. 
(KOO. 
(LM1. 
(MM  I. 
(NO.N 
C  TO  ( 

F(X(1) 
C  13 
IF(X( 

CONTI 

C  TO  1 

C  15 
IF(  X( 
CONTI 

C  TO  1 

C  17 
I F  (  X  ( 
IF(  Y( 
CONTI 

=KO+MO 

=  LC+1 

C  19 

K=K-y 
1  =  1-1 

U(KJ  = 
V<K)  = 

N  =  MO 

N=1.0/ 

CC-LOG 

5=MO+l 

C  80 

NES  TO 

TO 

I  F(  I  . 

X3  =  U( 

Y3  =  V( 

X4=U( 

Y4=V( 

/*3  =  X4 

E3=Y4 

IFtMC 

IF(LO 

£4=A3 

E4=B3 

GG  TO 

£2=A3 

/1  =  A2 

E2  =  B3 

E1=B2 

GC  TC 

>2=X3 

Y2=Y3 

X3  =  X4 


C-l 

1 

1 

LE.C) 

GE.3) 

LE.C) 

LE.C) 

E.LP1* 

litlfc) 


M0  + 
,  M 


1  ) 

DC 


GO 
GO 
GC 
GO 
GO 


TO 
TO 
TO 
TO 
TO 


9  0 
90 
91 
92 
93 


-X(2)  ) 

1  =  3, LO 

I-l)-X( 

NUE 

8 

1=3, LO 

I-l)-X( 

NUE 

8 

1  =  2, 

1-1  ) 

i-n 

NUE' 


I)  ) 


I)  ) 


12,95,14 
13,95,9  6 

96,95,15 


LO 
•  NE 

.EQ 


(  I)  ) 
(  I  )  ) 


GC 

GC 


TO  17 
TO  97 


J=1,L0 

0 

X{  I  ) 
Yd  ) 

RM 
P 

1  =  1, LO 
PICK  UP 
ESTIMATE 

GT.l) 

1) 

1) 

FO+1) 

M0.+  1J 

-X3 

-Y3 

Ml.EQ.C) 

.NE.2) 


(33,32)  , 
+A3-A4 
+A2-A3 
+B3-B4 
+B2-63 

(51,56) , 


NECESSARY  X  AND  Y  VALUES  MC 
TFEM  IF  NECESSARY 
GO  TO  40 


M3  = 
GO 


33/A3 
TO  41 


MCO 


MCO 
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Y3=Y4 

X4  =  X5 

Y4  =  Y5 

Ai=A2 

£1=B2 

A2  =  A3 

B2=B3 

A3  =  A4 

B3=B4 

I  F ( I .GE.LM1) 

GO 

TO  42 

41 

K5=K5+M0 

X5=U(K5) 
Y5=V( K5) 
A4=X5-X4 
E4=Y5-Y4 

IF(MCfl.EQ.O) 

M4  = 

=B4/A4 

GO  TO  43 

42 

IF(MDM1 .NE.O) 
B4=B3-i-83-B2 

A4  = 

=A3+A3-A2 

4^ 

IF(I  .EQ..1I 

GO 

TO  3  1 

GC  TO  (50,55) , 

MCO 

NUKE 

RICAL  DIFFERENTIATI 

ON 

50 

T2=T3 

51 

W2=DA8S(N4-M3) 

5t 


57 


58 


W3=0A6S(K2-M1 ) 

SW=W2+W3 

IF(SW  .KE.O.O) 

W2=0. 5 

W3=0.5 

SW=1.C 

T3=(k2*M2+k3*M3) 

IF(I-l)  30,00,60 

CCS2=C0S3 

SIN2=SIN3 

W2=CABS( A3-B4-A4 

V\3=DABS(  A1'*B2-A2 

IF(W2+W3.NE. 0.0) 

k2=DSCRT(A3*A3+B 

h3=DSQRT(A2*A2+6 

CCS3=h2*A2+W3*A3 

SIN3=W2*B2+W3*83 

R=COS3*C0S3+SIN3 

IF(R.  EC. 0.0) 

R=DSGRT(R) 

CCS3=CCS3/R 

SIN3=SIN3/R 


GC    TO    52 


/SW 


*B31 
*  B 1 ) 

GO    TO    5  7 
3*63) 
2*62) 


*  S I  N  3 

GO    TO    58 


IF( 1-1) 


C    DETERMINATION 
6C 


a  c ,  g  o ,  6  5 

OF    THE    C 


COEFFICIENTS 

iv2-T3)/A2 

)/(A2*A2) 


*B2) 

2+C0S2+C0S3) 


Q2=(2 .0*(M2~T2)+ 
Q3=(-M2-M2+T2+T3 

GC  TO  70 
65    R  =  DS0RT(A2J:A2-t-B2 
P1  =  R.*CGS2 
P2=3.C*A2-R#(C0S 
P3=A2-P1-P2 
C1=R*SIN2 

C2=3.C*B2-R*(SIN2+SIN2+SIN3) 
C3=B2-0l-Q2 
GC  TO  75 
COMPUTATION  OF  THE  POLYNOMIALS 
7C    DZ=A2*RM 
Z  =  0.0 

DO  71   J=1,MM1 
.  K  =  K  +  1 
Z=Z+DZ 
U(K)=PC+Z 
71     '  V(K)=C0+Z*(G1+Z*(G2+Z*Q3)) 

GO  TO  7S 
75    Z=0.0 

DO  76   J=1,MM1 
K  =  K  +  1 
Z=Z+RM 


119 


7 
7 

a 

82 


81 

C  EP 

9 

9 

9 


92 
95 
96 


9 
9 
9 

,  FC 
209 
209 
209 
209 
209 
2C9 
209 
209 

209 


6 

c 

• 

c 
c 

A 
N 
C 
A 
R 
RGB 
C  U 
G 

1  k 
G 

2  W 
G 

G 

£  k 

G 

6  h 
G 

7  W 

8  k 

9  M 
P 

RMA 
C  F 

1  F 

2  F 

3  F 

5  F 

6  F 

7  F 
g  F 


U(K 

V(K 

K  =  K+1 

CCNTI 
C  82  I 
NGL( I  ) 

M1=N0- 
C  81  K 
NGL(K) 
ETURN 
EXIT 
RITE  ( 
C  TO 
RITE 
C  TO 
RITE 
C  TO 
RITE 
C  TO 
RITE 
C  TO 
RITE 
0  TO 
RITE 
RITE 
RITE 
ETORN 
T  STAT 
CRMATf 
CRMATt 
CRMATt 
CRMAT  ( 
CRMATt 
CRMAT ( 
CPMAT( 
CRMATt 


)=PC+Z*tPl+Z*(P2+Z*P3)  ) 
)=QC+Z*(Q1+Z*(G2+Z*Q2)  ) 

NUE 

=  1,12 

=C.CDO 

1 

=12,NM1 

=DATAN2(  (UtK+l)-UtK-l)  I   , 


6,2C9G) 

9 

6,2C91) 

c 

6,2  092} 

9 

6,2C93) 

9 

6,2C95) 

8 

6,2  096) 

8 

6,2097) 

6T2C93) 

6,2C99) 


(  V(K  +  1) -VtK-1)  )  } 


I  ,Xt  I)  ,Y(  I) 

MDC,L0,MO,N0 


EVENTS 
1X/3LH 
1X/22H 
IX/22H 
LX/25H 
1X/27H 
1X/33H 
1X/33H 


1 

9  FCRMATt 
1 
2 
END 


7h 

7H 

3HM 
36H 


1 


*<*  ~t*  & 
:!«  $  if 

=  »I4 


MD  CUT  OF  PROPER 
L  =  1  OR  LESS./) 
M  =  1  OR  LESS./) 
IMPROPER  N  VALUE 


RANGE  ./) 


/) 


IDENTICAL  X  VALUES./) 


M  C  =  ,  1 4 


X  VALUES  OUT  OF 

IDENTICAL  X  AND 

1CX,6HX(I)  =.E12 

10X,6HY(I)  = , E 12 

8X,3HL  =,I5,8X, 


-,  I5,dX,3HN  =, 15/ 

ERROR  DETECTED  IN  ROUTINE 


SEGUED  E./) 
Y  VALUES./) 
3, 
3) 


CRVFIT) 
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CCRNER  NODAL  PCINTS  AND  THEIR  COORDINATES 


NCCE 

1 

X  = 

-3. £00 

R  = 

0.0 

NCDE 

3 

X  = 

-3. EGG 

R  = 

0.300 

NCCE 

K 

X  = 

-3. £00 

R  = 

C.600 

NCCE 

7 

x= 

-3  .800 

R  = 

0.800 

NCDE 

9 

x= 

-3. EGG 

R  = 

1.000 

NCDE 

IS 

x= 

-3  .300 

R  = 

0.0 

NCCE 

21 

x  = 

-3  .300 

R  = 

0.300 

NODE 

22 

x  = 

-3.2C0 

R  = 

0.600 

NCCE 

25 

x= 

-3.300 

R  = 

0.800 

NCDE 

27 

x  = 

-3.300 

R  = 

1.000 

NODE 

37 

x= 

-2.6C0 

R  = 

0.0 

NCCE 

39 

x  = 

-2  .800 

R  = 

0.300 

NODE 

41 

x= 

-2.600 

R  = 

0.6  00 

NODE 

43 

x= 

-2  .800 

R= 

C.800 

NCCE 

45 

x  = 

-2  .800 

R  = 

1  .000 

NODE 

55 

x  = 

-2.300 

R= 

0.0 

NODE 

57 

x= 

-2  .300 

R= 

0.300 

NCDE 

59 

x  = 

-2.20C 

R= 

0.600 

NCCE 

61 

x= 

-2.3CG 

R  = 

C.600 

NCCE 

63 

x= 

-2  .300 

R= 

1.000 

NCDE 

73 

x= 

-1  .9CC 

k  = 

0.0 

NCCE 

75 

x  = 

-1  .900 

R  = 

0.300 

NCCE 

77 

x  = 

-1  .9  00 

R= 

0.600 

MODE 

79 

x= 

-1.90C 

R  = 

G.800 

NCCE 

81 

x  = 

-1  .900 

R= 

1.000 

NODE 

91 

x  = 

-1.50C 

R  = 

0.0 

NCDE 

93 

x  = 

-1.500 

R  = 

C.3C0 

NCCE 

95 

x= 

-1  .500 

R  = 

0.600 

NCDE 

97 

x= 

-  1  .  5  0  C 

R  = 

C.8  00 

NCCE 

9  9 

x= 

-1.500 

R  = 

1.000 

NCCE 

110 

x= 

-1  .3  50 

R= 

1  .000 

NODE 

111 

x= 

-1.2C0 

R  = 

0.0 

NCCE 

113 

x  = 

-I  .200 

R= 

0.300 

NODE 

115 

x= 

-1  .200 

R  = 

0.600 

NODE 

117 

x= 

-1  .20C 

R  = 

0.800 

RCCE 

119 

x= 

-1  .200 

R  = 

0.97  0 

NODE 

121 

x  = 

-  1  . 2  C  C 

R  = 

1.000 

NODE 

125 

x  = 

-I  .  150 

R  = 

I. 000 

NCCE 

138 

x  = 

-1  .COO 

R  = 

0.0 

NODE 

14C 

x  = 

-l.CCG 

R  = 

0.300 

NCDE 

1-42 

x  = 

-1  .000 

R= 

G.60C 

NCDE 

144 

x  = 

-1  .COC 

R  = 

0.800 

NODE 

147 

x= 

-1.  100 

R  = 

C.970 

NCCE 

150 

x= 

-1  .100 

R  = 

1.000 

NCDE 

153 

x  = 

-1.C5C 

R  = 

1.000 

NCCE 

156 

x= 

-1  .100 

R  = 

0  .  9  7  0 

NCCE 

158 

x  = 

-1  .000 

R  = 

1.000 

NODE 

162 

x  = 

-C.950 

R  = 

I. COO 

NCDE 

173 

x  = 

-C  .ECC 

R= 

CO 

NCDE 

175 

x  = 

-C.60C 

R  = 

0.300 

NOCE 

177 

x  = 

-c.eco 

R= 

C.600 

NCCE 

179 

x= 

-0.800 

R  = 

0.800 

NODE 

1£2 

x  = 

-  C  .  9  C  0 

R  = 

0.955 

NOCE 

164 

x= 

-O.EOO 

R  = 

0.940 

KCCE 

189 

x  = 

-0.900 

R  = 

1.000 

NOCE 

19C 

x  = 

-C.E5G 

R  = 

1.000 

NCCE 

195 

x  = 

-0  .765 

R  = 

1.005 

NCDE 

2CC 

x  = 

-0.800 

R  = 

1  .000 

NCDE 

2C7 

x  = 

-C.6GC 

R  = 

O.C 

NCCE 

2C9 

x  = 

-C  .600 

R= 

0.300 

NCDE 

211 

x  = 

-C.6CC 

R  = 

0.600 

NCCE 

213 

x  = 

-0.630 

R  = 

0-814 
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NCDE 

217 

X  = 

-0  .725 

R  = 

0.952 

NODE 

219 

X  = 

-C.651 

R  = 

C.964 

NCCE 

224 

x= 

-0  .729 

R  = 

1.011 

NCDE 

226 

x  = 

-0.694 

R  = 

1  .019 

NCDE 

228 

x= 

-C.66C 

R  = 

1.028 

NCCE 

230 

x  = 

-0  .630 

R  = 

1  .037 

NCDE 

232 

x  = 

-C.601 

R  = 

1.046 

NCDE 

244 

x= 

-C  .400 

R  = 

0.0 

NCCE 

246 

x  = 

-0  .400 

R  = 

0.300 

NODE 

248 

x= 

-C.4CC 

R  = 

0.600 

NCCE 

250 

x= 

-0  .472 

R= 

0.h35 

NCDE 

253 

x= 

-0.587 

R= 

0.981 

NOCE 

255 

x= 

-C.522 

R  = 

C.999 

NCCE 

26  1 

x  = 

-C  .572 

R  = 

1.057 

NCDE 

263 

x  = 

-C.544 

R  = 

1  .069 

NODE 

26  5 

x= 

-C.521 

R= 

1.030 

NCDE 

267 

x  = 

-0  .498 

R  = 

1.092 

NODE 

279 

x  = 

C.C 

R  = 

0.0 

NCCE 

281 

x= 

-C  .  100 

R  = 

0.300 

NOCE 

233 

x= 

-0.160 

R  = 

O.bOO 

NCDE 

2S5 

x= 

-0.318 

R  = 

C.860 

NCCE 

288 

x  = 

-0  .468 

R  = 

1.020 

NODE 

29C 

x  = 

-C.414 

R  = 

1.042 

NCDE 

296 

x= 

-C.476 

R  = 

1.105 

NCCE 

298 

x  = 

-0.455 

R  = 

1.120 

NOCE 

3CC 

x  = 

-C.444 

R= 

1.128 

NCCE 

3C2 

x= 

-0  .433 

R  = 

1.13  7 

NOCE 

3C6 

x  = 

-0.378 

R  = 

1.059 

NODE 

3  16 

x= 

C.70C 

R  = 

0.0 

NCCE 

318 

x= 

C  .  4  5  C 

R= 

0  .270 

NCDE 

32C 

x  = 

C.25G 

R  = 

0.450 

NCCE 

3  22 

x= 

C  .  C  5  0 

R= 

C.620 

NCCE 

324 

x  = 

-0  .181 

R  = 

0.388 

NOCE 

3  28 

x  = 

-C.343 

R= 

1.076 

NCDE 

334 

x= 

-0  .422 

R  = 

1.146 

NODE 

336 

x  = 

-0.412 

R  = 

1.156 

NODE 

"3  °  6 

x= 

-C.404 

R  = 

1.  165 

NCCE 

340 

x  = 

-0  .395 

R= 

1.174 

NCDE 

357 

x  = 

1  .  15C 

R  = 

0.0 

MODE 

36  1 

x  = 

-0.387 

R  = 

1.183 

NCCE 

363 

X- 

-0  .360 

R= 

1  .193 

NCDE 

366 

x= 

-C.314 

R= 

1.104 

NODE 

36  9 

x= 

L  .  6  C  0 

R  = 

0.0 

NCCE 

371 

x  = 

1  .350 

R= 

0.200 

NODE 

373 

x  = 

1.4C0 

R  = 

0.350 

NCCE 

375 

x  = 

C  .74  0 

R= 

0.510 

NCCE 

377- 

x= 

0.50C 

R  = 

0  .c60 

NCCE 

379 

x= 

C.25C 

R  = 

0.790 

NCCE 

381 

x  = 

-0  .0  65 

R  = 

0.991 

NCDE 

383 

>  = 

-C.286 

R  = 

1.133 

NCCE 

386 

x  = 

-0.375 

R  = 

1.200 

NCCE 

388 

x  = 

-0  .371 

R  = 

1.206 

NOCE 

39  2 

x  = 

-C.275 

R  = 

1.155 

NCCE 

408 

x  = 

1  .59  0 

R  = 

O.L/9 

NODE 

412 

x  = 

-C  .367 

R  = 

1  .213 

NODE 

4  14 

x= 

-C.363 

R= 

1.22C 

NCCE 

420 

x= 

-0.C32 

R  = 

1.G75 

NCDE 

422 

x  = 

C.3CC 

R  = 

0.930 

NCDE 

423 

x= 

1.55C 

R  = 

0.397 

NCCE 

424 

x  = 

1  .470 

R  = 

0.632 

NOCE 

426 

x= 

1  .150 

R  = 

C.720 

NCCE 

426 

x= 

C  .£00 

R  = 

C.810 

NCCE 

43C 

x= 

0  .500 

R  = 

C.890 

NOCE 

432 

x= 

-0.26  4 

R  = 

1.177 

NCCE 

434 

x= 

-0  .358 

R  = 

1.229 

NCDE 

4  36 

x  = 

-C.354 

R  = 

1  .238 

NCDE 

4  5  5 

x= 

-C.256 

R  = 

1.205 

NCCE 

456 

x= 

-0  .023 

R  = 

1  .179 

NODE 

46C 

x  = 

-C.35C 

R= 

1.248 
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NCCE 

46  5 

X  = 

-0.24  6 

R=    1 

..257 

NCCE 

467 

x- 

-0  .249 

R=    1 

..233 

NODE 

4tS 

X  = 

C.200 

R  =    ] 

..100 

NCCE 

471 

X  = 

0  .550 

R=     1 

.100 

NCCE 

473 

x  = 

C.82C 

R=    1 

.100 

NODE 

475 

x= 

1  .  1 6  C 

R=    J 

l.  102 

NCCE 

4  77 

x  = 

1  .290 

R=    C 

).792 

NCDE 

476 

x  = 

L.28C 

R=    C 

:.960 

NCCE 

483 

x= 

I.C30 

R=    1 

L.224 

NCCE 

494 

x  = 

-0.250 

R=    ] 

..259 

NODE 

499 

x= 

-C.243 

R-    i 

..265 

NCCE 

501 

x  = 

-C  .240 

R= 

L.2  74 

NCDE 

5C3 

x^ 

C.675 

R=    i 

..340 

NOCE 

5C5 

x= 

-C.238 

R=    J 

1.282 

NCCE 

507 

x  = 

-0  .226 

R=    ] 

l.291 

NCDE 

512 

x= 

-C.C53 

R=     J 

..271 

NCCE 

5  14 

x= 

C.22C 

R=    J 

L  .250 

NCCE 

516 

x= 

0.480 

R=     ] 

..300 

NOCE 

518 

x= 

C.725 

R=    ] 

..426 

NCCE 

522 

x= 

-0.251 

R=    ] 

L  .285 

NCCE 

525 

x  = 

-0.335 

R=    ] 

..,'99 

NGDE 

529 

x  = 

-0.233 

R=    J 

L.307 

NCCE 

531 

x  = 

-0  .254 

R  = 

L.310 

NCDE 

54  1 

x  = 

C  .  5  9  C 

R=    ] 

..487 

NCCE 

54  2 

x= 

C  .  4  6 ■  C 

R  = 

l.532 

NCCE 

544 

x  = 

0.350 

R=    1 

L  .561 

NODE 

546 

x= 

C  .  1 6  C 

R=    ] 

..400 

NCCE 

54  8 

x= 

-0.C86 

R  = 

1.362 

NCCE 

553 

x  = 

-0.232 

R=    ] 

L  .316 

NODE 

C  C  £ 

x  = 

-0.231 

R  = 

L.324 

NCCE 

558 

x  = 

-0.257 

R  = 

L  .335 

NCDE 

561 

x  = 

C.24C 

R=    j 

L.582 

NCCE 

57C 

x= 

-C.264 

R  = 

L.352 

NCCE 

573 

x= 

-0.330 

R=    J 

L  .330 

NCDE 

5  76 

x= 

-C.23C 

R=    ] 

L.337 

NCCE 

57  9 

x= 

-0  .229 

R  = 

1.343 

NOCE 

585 

x= 

-0  .229 

R  = 

1.350 

NOCE 

587 

x= 

-C.271 

R  = 

L.37  0 

NCCE 

569 

x= 

-0  .  135 

R  = 

L.4  15 

NCDE 

591 

x  = 

C.C6C 

R=    i 

L.480 

NODE 

593 

x= 

C.  155 

R  = 

1.592 

NCCE 

595 

x= 

0.070 

R  = 

L.598 

NODE 

597 

x= 

c.c 

R=    ] 

L.600 

NCCE 

5  9  9 

x= 

0  .  C  2  5 

R  = 

L.60C 

NCCE 

lit 

x= 

-0.160 

R  = 

L.480 

NODE 

6C8 

x  = 

-C.212 

R  = 

1.414 

NCCE 

6  17 

x  = 

-C  .038 

R  = 

L  .640 

NCDE 

618 

x= 

-C.23C 

R=    ] 

L  .440 

NCCE 

62C 

x= 

-C.  177 

R  = 

L.591 

NCCE 

627 

X  = 

-0.C77 

K  = 

L  .679 

NODE 

626 

x= 

-C.26C 

R  = 

L  .520 

NCCE 

62  2 

x= 

-0  .249 

R  = 

L.675 

NCCE 

635 

x  = 

-0.C91 

R  = 

L.59  3 

NODE 

627 

x= 

-0.116 

R  = 

L.719 

NCCE 

628 

x  = 

-C  .4  10 

R  = 

L  .600 

NCDE 

64  2 

A  - 

-C.222 

R  = 

L  .758 

NCCE 

644 

x= 

-C.  166 

R  = 

1.675 

NCCE 

647 

x= 

-0.154 

R.= 

1.759 

NOCE 

646 

x  = 

-C.481 

R  = 

L.765 

NCCE 

654 

x  = 

-0  .295 

R  = 

1.84  1 

NOCE 

657 

x  = 

-0.231 

R  = 

L  .838 

NCCE 

658 

x= 

-C.551 

R  = 

1.851 

NCCE 

66  2 

x= 

-0  .468 

R  = 

1.923 

NCCE 

667 

x  = 

-C.308 

R  = 

L.V17 

NCCE 

666 

x= 

-C.621 

R  = 

1.936 

NCCE 

67C 

x  = 

-0  .542 

r=   ; 

I  .006 

NGDE 

672 

x  = 

-C.462 

r=   ; 

'.076 

NCCE 

677 

x  = 

-0  .285 

R  = 

1.996 
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